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Abstract. It is the purpose of this paper to rediscuss the 
circumstellar properties of S stars and to put these prop- 
erties in perspective with our current understanding of 
the evolutionary status of S stars, in particular the intrin- 
sic/extrinsic dichotomy. This dichotomy states that only 
Tc-rich ('intrinsic') S stars are genuine thermally-pulsing 
asymptotic giant branch stars, possibly involved in the M- 
S-C evolutionary sequence. Tc-poor S stars are referred to 
as 'extrinsic' S stars, because they are the cooler analogs 
of barium stars, and like them, owe their chemical pecu- 
liarities to mass transfer across their binary system. 

Accordingly, an extensive data set probing the circum- 
stellar environment of S stars (IRAS flux densities, maser 
emission, CO rotational lines) has been collected and crit- 
ically evaluated. This data set combines new observations 
(9 stars have been observed in the CO J = 2 — 1 line and 
3 in the CO J = 3 — 2 line, with four new detections) 
with existing material (all CO and maser observations of 
S stars published in the literature). The IRAS flux den- 
sities of S stars have been re-evaluated by co-adding the 
individual scans, in order to better handle the intrinsic 
variability of these stars in the IRAS bands, and possible 
contamination by Galactic cirrus. 

In the (if — [12], [25] — [60]) color-color diagram, S stars 
segregate into five distinct regions according to their Tc 
content and ZrO/TiO, C/0 and IR spectral indices. Stars 
with photospheric colors (populating 'Region A') may be 
identified with extrinsic S stars. For the other regions char- 
acterized by different excess levels in the 12, 25 and 60 fim 
bands, several diagnostics (like the IRAS spectral class, 
maser emission, and shape of CO rotational lines) have 
been collected to infer the physical properties of the dust 
shell. A simple radiative-transfer code has also been used 
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to infer the chemical nature (carbonaceous or silicate) of 
the dust grains from the observed IR colors. S stars with 
large K — [12] excesses and moderate [25] — [60] excesses 
(populating Regions B and C) exhibit the signatures of 
oxygen- rich shells (9.7 fim silicate emission and SiO maser 
emission). The situation is less clear for S stars with small 
K — [12] and moderate [25] — [60] indices (populating Re- 
gions D and E). Their IR colors are consistent with car- 
bonaceous grains (as is their featureless IRAS spectrum, 
and absence of silicate or SiO maser emission) , but these 
features may equally well be explained by a detached shell. 
For many of these stars with a large 60 /im excess, the shell 
is indeed resolved by the IRAS beam at 60 /im. The pro- 
totypical SC star FU Mon is among these. Since SC stars 
are believed to be in a very short-lived evolutionary phase 
where C/0 = 1 within 1%, FU Mon may be a good candi- 
date for the 'interrupted mass-loss' scenario advocated by 
Willems & de Jong (1988). The CO hue profile of FU Mon 
is also peculiar in being quite narrow (Ve — 2.8 km s^^), 
suggesting that the mass loss has just resumed in this star. 

Mass loss rates or upper limits have been derived for all 
S stars observed in the CO rotational lines, and range from 
< 2 10"^ Mq y-i for extrinsic S stars to 1 10"^ Mq y-^ 
(the Mira S star W Aql). These mass-loss rates correlate 
well with the K — [12] color index, which probes the dust 
loss rate, provided that M > 10~^ Mq y~^. Small mass- 
loss rates are found for extrinsic S stars, consistent with 
their not being so evolved (RGB or Early-AGB) as the 
Tc-rich S stars. This result does not support the claim 
often made in relation with symbiotic stars that binarity 
strongly enhances the mass-loss rate. 

Key vifords: Stars: mass-loss - Stars: AGB - Stars: S - 
Stars: late-type - Infrared: stars - Radio lines: stars 
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1. Introduction 

The S stars arc latc-typc giants whose spectra resem- 
ble those of M giants, with the addition of distinctive 
molecular bands of ZrO (Merrill 1922). The presence of 
ZrO bands is often considered as a direct consequence of 
molecular equilibrium in the special circumstances where 
the atmospheric C/0 ratio is within 10% of unity (e.g. 
Scalo & Ross 1976). However, Piccirillo (1980) has shown 
that the above statement is only valid for stars with 
T < 3000 K. At higher temperatures, an enhanced Zr 
abimdance, rather than a C/0 ratio close to unity, is 
the dominant factor in the development of strong ZrO 
bands. Detailed abundance analyses (Smith & Lambert 
1990) have shown that the envelopes of S stars are en- 
riched in heavy elements like Zr, and so bear the signature 
of the s-process of nucleosynthesis (Kappeler et al. 1989). 
Although larger than in M giants, the C/0 ratio of S stars 
is not necessarily close to unity (Smith & Lambert 1990), 
except in the so-called SC stars (Dominy et al. 1986). 

When S stars were still believed to be objects with 
C/0 close to unity, they were naturally considered as tran- 
sition objects between M giants (C/0 < 1) and carbon 
stars (C/0 > 1) on the asymptotic giant branch (AGB) 
(Iben & Renzini 1983). Support to this scenario is pro- 
vided by observations of S stars on the upper AGB of 
globular clusters in the Magellanic Clouds (Bessell et al. 
1983; Lloyd Evans 1983, 1984, 1985). In this evolution- 
ary phase, low- and intermediate-mass stars are charac- 
terized by a double (H, He) burning-shell structure which 
is thermally unstable. The thermal instabilities ('thermal 
pulses') developing in the He-burning shell are the site 
of a rich nucleosynthesis (Frost & Lattanzio 1995), prob- 
ably including the s-process, although its detailed mode 
of operation remains poorly understood (e.g. Sackmann 
& Boothroyd 1991; Herwig et al. 1997). In the receding 
phase of the thermal instability, the convective outer en- 
velope may plunge ('third dredge-up') into the intershell 
zone containing the He- burning ashes, and bring fresh car- 
bon and s-process elements to the surface. 

However, several observations have challenged this tra- 
ditional M S C evolution sequence. The first set of ob- 
servations relates to Tc, an element with no stable iso- 
topes, discovered in the spectra of some S stars by Mer- 
rill (1952). If the s-process indeed occurred during recent 
thermal pulses in S stars, Tc should be observed at the 
surface along with the other s-process elements (Mathews 
et al. 1986). Little et al. (1987) found however that only 
long-period Mira or semiregular S stars (i.e., intrinsically 
bright S stars) exhibit Tc lines. Second, the broad range 
of IRAS colors exhibited by S stars (e.g., Jorissen et al. 
1993) is difficult to reconcile with the idea that they rep- 
resent a brief transition phase as the star evolves from 
an oxygen- rich M giant with C/0 < 1 into a C star with 
C/0 > 1. M and C stars occupy well-defined regions in the 
IRAS color-color diagram, and it is not clear how S stars 



fit into the (much debated) evolutionary sequence joining 
M stars to C stars in that diagram (Willems & de Jong 
1986, 1988; Chan & Kwok 1988; Zuckerman & Maddalena 
1989; de Jong 1989). 

These problems received a new impetus with the dis- 
covery that the barium stars, a family of peculiar red giant 
(PRC) stars of spectral type G and K, are all members 
of binary systems (McClure et al. 1980; McClure 1983). 
Iben & Renzini (1983) were the first to propose that Tc- 
poor S stars could perhaps be the cooler analogs of the 
barium stars. Long-term radial-velocity monitorings con- 
firmed this suggestion, and it is now clear that Tc-poor S 
stars are binary stars (Smith & Lambert 1988; Brown et 
al. 1990; Jorissen & Mayor 1992; Jorissen et al. 1993; John- 
son et al. 1993) with orbital elements identical to those of 
barium stars (Jorissen et al. 1997). Tc-poor S stars are 
now referred to as 'extrinsic S stars', because, like barium 
stars, they owe their chemical peculiarities to mass trans- 
fer across the binary system. On the contrary, Tc-rich, 'in- 
trinsic S stars' are genuine thermally-pulsing (TP) stars 
on the TP- AGB. Since the C/0 ratio of extrinsic S stars 
depends on the details of the mass accretion process, it 
is not necessarily close to unity, but as discussed above, 
neither abundance analyses nor predictions of molecular 
chemical equilibrium really require C/0 to be close to 
unity in S stars. 

In Paper I (Jorissen et al. 1993: sec; also Groenewegcn 
1993), it was shown that the correlation Tc-poor/binary 
found for S stars could be extended to their IRAS col- 
ors, since another distinctive property of binary, Tc-poor 
S stars is the absence of IR excesses. Their IRAS colors 
simply reflect the photospheric colors, contrary to Tc-rich 
S stars which usually exhibit IR excesses. These IR ex- 
cesses are caused by circumstellar dust, and are indicative 
of substantial mass loss, thus suggesting that Tc-rich S 
stars are more massive and/or more evolved than Tc-poor 
S stars. 

The possibility offered by the circumstellar properties 
of S stars to probe their evolutionary status has been 
used in several recent studies (Jura 1988; Chen & Kwok 
1993; Bieging & Latter 1994; Sahai fc Liechti 1995). How- 
ever, these studies still rc;ly on hypotheses not fully consis- 
tent with the dichotomy recently found among S stars, as 
discussed above, and their conclusions may therefore be 
somewhat biased. For example, it makes no sense to test 
the AGB evolutionary sequence M-S-C using a sample of 
S stars not properly cleaned from its extrinsic content. In 
addition, models inferring the chemical nature (carbona- 
ceous or silicate) of the dust grains from models assuming 
that C/0 is close to unity in the photosphere do not sam- 
ple the whole parameter space occupied by these stars. 
A similarly incorrect corollary consists of concluding that 
dust production, and thus mass loss, is not very efficient 
in S stars because dust-seed molecules no longer form in 
the absence of any free C or O atoms, since these are all 
tied up in CO when C/0 ~ 1. 
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It is the purpose of this paper to rediscuss the cir- 
cumstellar properties of S stars and to put these prop- 
erties in perspective with our current understanding of 
the evolutionary status of S stars, in particular the intrin- 
sic/extrinsic dichotomy. An extensive data set probing the 
circumstellar environment of S stars (IRAS colors, maser 
emission, CO rotational lines) has therefore been collected 
and critically evaluated. This data set combines new ma- 
terial with existing results collected from the literature. 

The IRAS colors of S stars provide a first way to probe 
their circumstellar environment, and more specifically, to 
evaluate the amount and nature of the dust surround- 
ing these objects. A re-evaluation of the flux densities 
listed in the IRAS Point Source Catalogue (IRAS Science 
Team 1988; PSC) was necessary to circumvent the prob- 
lems of interpretation related to the intrinsic IR variabil- 
ity of these sources and to their possible contamination 
by Galac'tic cirrus emission. These effects are not always 
properly handled in the PSC. This re-evaluation was per- 
formed by co-adding the raw scans (Sect. 2.2). The clean 
flux densities were then used to define five regions in the 
{K - [12], [25] - [60]) color-color diagram (Sect. 2.3) which 
contain stars of similar extrinsic or intrinsic nature and of 
similar ZrO/TiO, C/0 and IR spectral indices. Several S 
stars with shells resolved at 60 (and sometimes 100) fxm 
have been found by comparing the source profile with the 
IRAS point source response function (Sect. 2.4). Resolved 
shells at 60 /xm appear to be correlated with large 60 fxm 
excesses, suggesting that these shells have detached from 
their parent star. A simple model of the dust shell has 
been used to predict its IR colors and to infer the chem- 
ical nature of the dust grains (Sect. 4), using constraints 
provided by the detection or non-detection of maser emis- 
sion (Sect. 3). Finally, mass loss rates have been derived 
in an homogeneous way (Sect. 5.3) from the intensities 
of the CO millimeter-wave lines, derived from new ob- 
servations with the Caltech Submillimeter Observatory 
(Sect. 5.1) or collected from the literature (Sect. 5.2). The 
mass loss rates, wind expansion velocities, IR colors and 
extrinsic/intrinsic nature are then discussed together in 
Sect. 6. 



2. Infrared colors of S stars 
2.1. The sample 

The sample of S stars considered in this paper was se- 
lected from the list of Chen et al. (1995), which provides 
associations between IRAS sources from the PSC and 
S stars from the General Catalogue of Galactic S Stars 
(GCGSS, Stephenson 1984) and from an additional later 
list (Stephenson 1990). Among these, only stars having 
flux densities of good quality (i.e. with a quality flag of 
3 in the PSC) in the 12, 25 and 60 /xm bands have been 
retained. 



Several stars that have probably been misclassified as 

5 stars were removed from the final sample. For exam- 
ple, stars from the original Westerlund & Olandcr (1978) 
sample were later recognized by Lloyd Evans & Catch- 
pole (1989) as actually being heavily-reddened M giants 
or supergiants. A few other cases of M giants or super- 
giants misclassified as S stars were identified by Keenan 

6 McNeil (1989; HR 3296 = GCGSS 500) and Winfrey et 
al. (1994; GCGSS 1314 and star 41 in Stephenson 1990). 
A detailed heavy-element abundance analysis (Lambert 
et al. 1995) has shown that the stars DE Leo (HR 4088) 
and HR 7442. although often considered as S stars, have 
normal abundances. According to Meadows et al. (1987), 
GCGSS 886 (IRAS 15194-5115) is now classified as a car- 
bon star and thus has not been retained in the final sample 
of S stars. Our final sample may still be somewhat con- 
taminated by M supergiants misclassified as S stars and by 
a few carbon stars. The star T Get for example was clas- 
sified as M5-6Se in the original paper by Keenan (1954) 
defining the S class, but it was reclassified as M5/M6Ib/II 
in the Michigan Spectral Survey (Houk & Cowley 1975). 
Since at the low plate dispersions used in classification 
work, these two types of spectra look similar (e.g. Lloyd 
Evans & Catchpole 1989), T Get has been kept in our 
final list until higher resolution spectra resolve these con- 
flicting classifications. The same holds true for TV Dra 
and OP Her (sec Table 2 of GCGSS). Note, however, that 
SC stars like RZ Peg, FU Mon and GP Ori are retained in 
the sample, since they may provide important clues to the 
evolutionary status of S stars as a whole. Similarly inter- 
esting are the two CS stars TT Gen and BH Cru that are 
known to exhibit ZrO bands at some times and C2 bands 
at others (Stephenson 1973; Lloyd Evans 1985). 

2.2. The IRAS flux densities 

The final sample consists of 124 S stars having fiux densi- 
ties at 12, 25 and 60 ^m flagged as being of good quality in 
the PSC. These stars are listed in Table 1. However, PSC 
flux densities suffer from several shortcomings that make 
them inadequate for the present study. First, the 60 /zm 
flux density may in some cases be seriously contaminated 
by Galactic cirrus emission, as shown by Ivezic & Elitzur 
(1995). Second, the PSC flux densities are not appropri- 
ate in the case of slightly extended or variable sources, as 
several S stars appear to be. Moreover, hysteresis of the 
detectors hampers the search for extended sources and 
should be properly identified. Finally, the detectors may 
saturate on very bright sources, those with flux densities 
in excess of 1000 Jy (like x Cyg, see Appendix A). 

In order to correctly handle these effects, the raw 
IRAS data for all S stars in our sample were reprocessed 
through the ADDSCAN procedure provided by the In- 
frared Processing and Analysis Center (IPAC^). This pro- 

^ IPAC is operated by the Jet Propulsion Laboratory and Cal- 
ifornia Institute of Technology for NASA 
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cedure has been used to co-add all scans passing within 
1'.7 of the target position. Before being co-added, the 
raw data are first interpolated using cubic splines and re- 
sampled at 10 data points per arcminute in all bands. 
A baseline is then defined for each individual scan by 
fitting a parabola to the data in a window extending 
from Tin to rout from the target (nn = ±2', ±2', ±2f5, ±4' 
and rout = ±7.'5, ±7'.5, ±10', ±13f5 in the 12, 25, 60 and 
100 /im bands, respectively). The r.m.s. residual ai of the 
data around the baseline fit is an indication of the back- 
ground noise in a given scan i (note that possible nearby 
sources with a peak flux density exceeding 2.5(Tj arc au- 
tomatically removed and do not enter the final residual 
calculation). The different scans are then co-added with 
weighting factors equal to l/trf. The noise ct in the co- 
added scan, computed in a similar way as for the indi- 
vidual scans, is smaller than that in any individual scan, 
and the reduction may be substantial when scans with 
very different orientations, sampling different regions of 
the sky around the target, are available, since structure 
in the Galactic cirrus emission contributes to the noise. 
Finally, a point-source template with adjustable width is 
fitted to the co-added scan in the target window irin. 
Three different estimates of the flux densities in each band 
are then computed: the peak flux density Fp, the template 
flux density and the 'zero- crossing' flux density F^. The 
template flux density Ft is derived from the template fit to 
the data, and generally agrees with Fp, unless the source 
is extended or hysteresis eff'ects are important. For ex- 
tended sources, the 'zero-crossing' flux density F,^ must 
be preferred. It corresponds to the integrated flux density 
between the 'zero crossings', which are defined as the first 
locations, moving outwards from the peak location, where 
the source profile intersects the baseline. In Sect. 2.4, a 
criterion based on the comparison of F^ and Fp has been 
designed to identify sources with possibly resolved shells. 
If that criterion is met, the F^ fiux density (identified by 
a '-I-' in Table 1) has been adopted instead of Ft. In the 
case of bright sources, detector hysteresis results in a trail 
extending along the scan direction, thus disturbing the 
template fit. In that case, Fp has been adopted. The tem- 
plate fit to sources embedded in strong Galactic cirrus 
emission generally resulted in an abnormally narrow pro- 
file, because the adopted baseline is then too high with 
respect to the base of the source signal. These contami- 
nated sources were thus readily identified by their narrow 
template fit, and were rejected. Finally, individual scans 
clearly contaminated by nearby sources have been elimi- 
nated. 

The adopted IRAS flux densities are listed in Table 1 

in columns 5 8. The 2.2 /im flux density is listed in col- 
umn 3. The calibration of Beckwith et al. (1976; 620 Jy 
corresponds to K = 0) has been used to derive the 2.2 /itm 
flux density from the K magnitude provided by different 
authors, as listed in column 4: 1: Catchpole et al. (1979); 
2: Neugebauer & Leighton (1969; Two-Micron Sky Sur- 
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Fig. 1. Comparison between the PSC flux densities and those 

derived from the ADDSCAN procedure (sec text) in the 12, 
25 and 60 /xm bands. Flux densities of variable sources at two 
different epochs are connected by a vertical line. At 60 fim, 
open circles correspond to extended sources (see Sect. 2.4) 



vey); 3: Mendoza & Johnson (1965); 4: Price (1968); 5: 
Chen et al. (1988); 6: Noguchi et al. (1991); 7: Gughelmo 
et al. (1993); 8: Epchtein et al. (1990); 9: Epchtein et al. 
(1987); 10: Lloyd Evans & Catchpole (1989). 

The Tc content (from Paper I, and Lambert et al. 
1995) is listed in column 9. The columns labelled 'LRS' 
and 'VC provide the classiflcation of the IRAS low- 
resolution spectrum, according to the IRAS Low Resolu- 
tion Spectrometer Catalogue (Olnon et al. 1986) or to the 
Volk & Cohen (1989) schemes, respectively. The optical 
spectral type is from the GCGSS. The variability type, 
the period and the range of variation of the magnitude 
(lowest minimum - highest maximum) are from the Gen- 
eral Catalogue of Variable Stars (Kholopov et al. 1985; 
GCVS). The variable name and/or the HD/HR designa- 
tions, when available, are given in the last column. If maser 
emission has been detected (see Sect. 3 and Table 3), the 
maser type (SiO or OH) is given between parentheses after 
the stellar designation. 

Several sources observed by IRAS a few months apart 
turn out to be strongly variable in the IRAS bands. For 
these, the ADDSCAN procedure has been run separately 
on the two groups of data, and the flux densities are listed 
on two separate lines in Table 1. The corresponding ap- 
proximate Julian dates have been derived from the 'Satel- 
lite Operation Plan' number attached to the scans and 
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from the mission chronology provided in Table III.C.l of 
the IRAS Explanatory Supplement (1988). 

The flux densities derived from the ADDSCAN pro- 
cedure generally agree with the PSC flux densities within 
20%, as shown in Fig. 1. In all bands, the scatter is larger 
at low flux densities. Several stars show much larger de- 
viations because of intrinsic variability. For these, it turns 
out that the PSC flux densities correspond to one epoch, 
whereas the flux densities at the other epoch may differ 
by as much as a factor of 2. An example of a strongly 
variable source, % Cyg, is discussed in Appendix A. In 
the 60 /im band, large discrepancies are also found for ex- 
tended sources, as expected (see Sect. 2.4). 

2.3. The infrared color-color diagrams 

The use of the ([12] - [25], [25] - [60]) diagram to probe 
the circumstellar shells surrounding late-type stars was 
first demonstrated by Hacking et al. (1985) and by van 
der Veen & Habing (1988, VH). VH defined defined var- 
ious regions (labelled I to VII) corresponding to circum- 
stellar shells with relatively homogeneous properties. Fig- 
ure 2 presents the ([12] - [25], [25] - [60]) diagram for 
the sample of S stars defined in Sect. 2.1. In this pa- 
per, the color index [i] — [j] is defined as [i] — [j] = 
-2.5log[{F{i)/Fo{i))/{F{j)/Fo{j))] (where F{i) refers to 
the non-color corrected flux density in the IRAS band i, 
and i^o(«) is a normalization flux density as given in the 
IRAS Explanatory Supplement 1988). With this normal- 
ization, black bodies in the Rayleigh- Jeans limit have a 
color index of 0. The colors for stars whose IRAS flux 
densities show large variations are plotted as two symbols 
joined by a line segment. 

There is a strong correlation between the position of 
an S star in the color-color diagram and the intensity of 
the spectral features distinctive of the S star class (as mea- 
sured by the ZrO/TiO or C/0 spectral indices in the spec- 
tral classification schemes devised by Keenan 1954 and 
Keenan & Boeshaar 1980, respectively: both will be called 
'ZrO index' for simplicity in the following). In particular, 
S stars with ZrO indices larger than 3 are mainly found in 
Regions VII and Via with a few more in Region III, while S 
stars with ZrO indices smaller than or equal to 2 are found 
mainly in Regions I and II. This segregation of S stars in 
the ([12] — [25], [25] — [60]) diagram is a clear indication 
of the inhomogeneous nature of this family of peculiar red 
giants. In the center of Region I lie the S stars with pho- 
tospheric colors. These arc in fact the binary and Tc-poor 
S stars ('extrinsic' S stars) that owe their chemical pe- 
culiarities to mass transfer in a binary system (Sect. 1). 
Another group of weak S stars is found at the boundary 
between Regions I and II, a zone generally devoid of stars 
in the ([12] - [25], [25] - [60]) diagram (see e.g. VH, and 
Lewis 1989). This group includes the prototypical Tc-rich 
Mira S variable x Cyg. The few S stars located in Region 
II, defined by VH as comprising 0-rich stars with 'young' 




0.2 0.4 0.6 0.8 
[12]-[25] 

Fig. 2. The ([12] - [25], [25] - [60]) diagram for S stars. The 
size of the square is proportional to the ZrO/TiO or C/0 spec- 
tral index of the S star (whichever is available), in the spectral 
classification scheme of Keenan (1954) or Keenan & Boeshaar 
(1980), respectively. MS and SC stars were assigned classes 1 
and 9, respectively. Crosses correspond to stars with no spec- 
tral type available in the GCGSS. Tc-poor S stars are repre- 
sented by filled squares, and colors of variable stars measured 
at different epochs are connected by a line segment 



circ\.mistellar shells, indeed have small ZrO indices. Most 
S stars populate Region VII, where C-rich circumstellar 
shells are found according to VH, and it is therefore not 
surprising that those S stars have large ZrO indices. Fi- 
nally, several stars (among which are many of type SC) 
are located in Regions Via and b, many of them having 
conspicuously resolved circumstellar shells (sec Sect. 2.4). 

Several authors (see e.g. Habing 1996) have argued 
that the K — [12] index is superior to the [12] — [25] in- 
dex for tracing the mass loss rate, the focus of the present 
study, because of its greater wavelength range and the fact 
that the photosphere emits more strongly at K while the 
shell emits more strongly at 12 fxm. Therefore, we felt that 
it was more meaningful to use the {K — [12], [25] — [60]) 
diagram to define groups of S stars with homogeneous IR 
properties, as follows (Fig. 3): 

Region A: S stars with photospheric IR colors ('extrinsic 
S stars'); 

Region B: S stars with no excesses at 60 /im, and small 
ZrO indices; 

Region C: S stars with excesses in all three 12 /xm, 25 /xm 
and 60 /im bands, and large ZrO indices; 
Region D: S stars with large ZrO indices and 60 /xm ex- 
cesses, but small 12 /im excesses; 
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Fig. 3. Same as Fig. 2 for the {K - [12], [25] - [60]) diagram. 
Note that the K — [12] index may be inaccurate for variable 
stars, since their K and [12] magnitudes were not measured 
simultaneously 



Region E: mainly SC stars with large 60 /um excesses and 
often resolved shells. 

The stars classified this way are plotted in the ([12] 
— [25], [25] — [60]) color-color diagram in Fig. 4, and it 
is seen that the above regions almost exactly correspond 
to those of VH, with the exception of Region B which 
encompasses both Regions I and II, and Region VII which 
is a blend of Regions C and D (though stars from Region 
D are concentrated in the upper left of Region VII). The 
motivation for creating Region D is apparent in Fig. 5, 
which presents the distribution of the various IR spectral 
types (as defined by Volk & Cohen 1989) across the {K — 
[12], [25] — [60]) diagram. Region D differs from Region 
C in having many stars exhibiting a stellar continuum in 
the IR (class S) and none with silicate emission (class E). 
As it will be shown in Sect. 4, these features suggest that 
the circumstellar shells in Region D contain C-rich grains, 
whereas silicate grains are found in the circumstellar shells 
of Region C. Note that Regions C and D are also quite 
distinct in the (K - [12], [12] - [100]) diagram (Fig. 6). 

Finally, let us remark that the different regions defined 
above have been denoted by uppercase letters to avoid 
confusion with the Regions a to d defined in Paper I from 
the lower-accuracy PSC flux densities. The new classifica- 
tion into Regions A-E is actually not very different from 
that in Paper I (matching uppercase with lowercase let- 
ters, and with Region d splitting into D and E), but it 
is clearly superior since several stars that appeared excep- 



Fig. 4. Same as Fig. 2, where the stars have been drawn with 
a symbol corresponding to the region they belong to in the 

(A'— [12], [25] — [60]) diagram (Fig. 3), as follows: open squares: 
Region A; filled circles: Region B; filled triangles: Region C; 
crosses: Region D; open circles: Region E 




Fig. 5. Same as Fig. 3, with symbols referring to the spectral 
class of the IRAS low-resolution spectrum as defined by Volk 
& Cohen (1989): S: stellar continuum; F: featureless spectrum; 
E: silicate emission; C: SiC emission; ?: not available. Spectral 
type assignments for S stars are from Chen et al. (1995) 
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Fig. 6. The {K - [12], [12] - [100]) diagram. Symbols are as in 
Fig. 4 



tional in Paper I fit well into the present classification. For 
instance, x Cyg, NQ Pup and Ori were exceptional as 
being Tc-rich stars in Region a. With the more accurate 
ADDSCAN flux densities, x Cyg moves to Region B (see 
also Appendix A) while NQ Pup and Ori move to Re- 
gion E. The only remaining outliers in this respect are HR 
Peg, the only Tc-rich star in Region A, and DY Gem, the 
only Tc-poor star not in Region A (but in C). DY Gem 
is however exceptional in many other respects: it is a very 
cool star (S8,5 corresponding to Toff = 3000 K; Smith & 
Lambert 1990) and a SRa variable with a very long pe- 
riod (1145 d). Moreover, it has the largest [12] — [100] index 
among stars in Region C (see Fig. 6). 

2.4- S stars with envelopes resolved at 60 fjm 

The possibility that some S stars may have circumstellar 
shells resolved by IRAS is now examined. As argued by 
Young et al. (1993a; YPK), the 60 /im band is best suited 
for that purpose, because it is not contaminated by Galac- 
tic cirrus emission as severely as is the 100 /xm band. The 
60 fim co-added scans were examined, and characterized 
as follows. First, the width of the template profile fitted to 
the source is compared to that expected for point sources, 
namely 2'.05 and 1.'44 full widths at the 25% and 50% 
levels, respectively (Levine et al. 1993). However, this cri- 
terion is not sensitive to resolved shells showing up as a 
weak plateau at the base of the profile. A simple crite- 
rion has therefore been devised, based on the compari- 
son of the 'peak flux density' Fp with the 'zero-crossing 



flux density' i^z (see Sect. 2.2). These quantities are stan- 
dard outputs from the ADDSCAN procedure. In the case 
of a point source, Fz and Fp arc identical within a few 
times the noise, measured as the rms a of the residuals 
along the baseline outside the signal range (i.e. between 
2.'5 and 10' from the target position in either directions) 
after baseline subtraction (see Sect. 2.2). In the case of 
an extended source, the fraction e of flux density in ex- 
cess of that of a point source is expressed by e = (F^ — 
Fp)/Fp. For estimating the significance of this excess, one 
has to be aware of the following effects. First, very bright 
(> 500 Jy) sources may have a characteristic six-pointed 
star shape due to reflection from the telescope secondary 
mirror struts. Since approximately 5% of the peak flux 
density may be contained in the star pattern, Fz/ Fp values 
of the order of 1.05 may be of instrumental origin (Levine 
et al. 1993). Bright sources are also affected by hysteresis 
in the detector which causes a trail in the signal along the 
scan direction in the outgoing part of the scan. This trail is 
easily recognized on individual scans as it causes an asym- 
metry in the template fit. However, if the co-added scan 
results from individual scans made in opposite directions, 
as is often the case, trails will be present in both directions 
and will mimick an extended plateau at the base of the 
profile. For bright sources (> 100 Jy), this effect may thus 
also cause spurious excesses of the order of a few percent. 
Finally, for fainter sources, an inhomogeneous background 
may also cause spurious detections. In this case, a way 
to estimate the significance of the excess e is to compare 
it with the inverse signal-to-noise ratio 1/SNR = rr/Fp, 
which is nothing more than the relative flux density excess 
expected from the background noise. The significance, or 
'quality factor' QF, of the observed flux density excess 
can then be expressed as QF = e Fp/a = e SNR, so that 
QF > 5 for a detection at the 5a level. Figure 7 presents 
the ratio F^/Fp vs SNR for S stars not confused by nearby 
sources, and several stars with QF > 5 are found. One has 
to be aware, however, that the co-addition process will be 
mostly effective in lowering the baseline noise far away 
from the target, where scans with different orientations 
sample different regions of the sky. But since they all in- 
tersect approximately on the target, the noise-averaging 
process will not be as effective on the target. A criterion 
based on the absolute noise level present on-target has 
therefore been considered as well, by requiring that, to be 
considered significant, an excess Fz — Fp should not be 
smaller than some given threshold value of the order of the 
background fluctuations in the 60 /im band. It was found 
that meaningful results are obtained by adopting 0.3 Jy 
as the typical background fluctuation on-target, combined 
with a quality factor of 5. In Fig. 7, sources satisfying the 
T'z — i^p > 0.3 Jy criterion have been represented by open 
squares. Since very few square symbols are located below 
the dashed line in Fig. 7, both criteria are fulfilled simulta- 
neously for most of the stars, and thus provide consistent 
conclusions regarding the resolved nature of these sources. 
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Fig. 7. The ratio F^/Fp between the 'zero-crossing flux den- 
sity' and the 'peak flux density' in the 60 /im band, vs the S/N 
ratio along the baseline (see text for details). Stars lying above 
the dashed line have a flux density excess in the 60 /xm band 
with a quality factor larger than 5. Stars with a flux density 
excess above the 0.3 Jy threshold are represented by squares. 
Squares lying above the dashed line thus satisfy both criteria 
defined in the text and are probably truly resolved sources. 
Filled squares and crosses correspond to sources flagged by 
YPK as extended or non-extended, respectively, at 60 /xm. Ob- 
servations of the same star at two different epochs are con- 
nected by a line segment. Sources confused with a nearby 
source are not plotted 

Sources flagged as extended at 60 /xm by the above 
criteria are listed in Table 2. Columns 1, 2 and 3 give the 
IRAS name, the variable star name and variability typo 
when available, respectively. The zero-crossing flux density 

is in column 4. Column 5 lists the flux density ratio 
Fz/Fp, and columns 6 and 7 the full widths W2b and 
at the 25% and 50% levels, respectively, to be compared 
with 2^05 and 1^44 for a point source (Levine et al. 1993). 
The quality factor QF is listed in column 8. Mass loss 
rates and wind terminal velocities are given in columns 9 
and 10 (see Sect. 5). 

A similar search for latc-typc giants with resolved 
shells was performed by Young et al. (1993a,b). These 
authors used a more sophisticated method based on the 
possibility of successfully fitting the signal by a point 
source surrounded by a circumstellar shell having 'reason- 
able' properties. Our simpler approach has the advantage 
of being applicable to fainter stars, and several have in- 
deed been detected, as seen in Fig. 8. However, for bright 
sources, our method is more vulnerable to spurious de- 



tections due to hysteresis (YPK used only the scan data 
taken prior to passing over the source). 

As can be scon in Table 2 (in column 11, YPK+ and 
YPK— denote sources flagged by YPK as extended or not 
at 60 yum, respectively), the two methods give conflicting 
results for four stars, R And, S Cas, W Aql and T Get 
among the 11 bright objects common to the two sam- 
ples. In the first three cases, our detection is probably 
an artefact due to detector hysteresis, since individual (as 
opposed to co-added) scans show an extended tail only 
on the side posterior to the passage over the source (see 
Fig. 10 below). The situation is less clear for T Get, as its 
100 /Ltm profile is wider than the point-source template, 
suggesting that this source may be truly extended. 

Figure 8 presents the flux density excess F.^/ Fp vs 
F(2.2), the flux density at 2.2 /xm, and reveals that the 
properties of the resolved shells in Regions B and C are 
very different from those of Regions D and E. In Regions B 
and C, the flux density excess is of the order of a few per- 
cent, with a maximum of 15% for Y Lyn. Because the flux 
density excess is so small, the envelopes around stars in 
Regions B and G can be resolved only for the stars closest 
to the sun [i.e. with large F(2.2)], as shown by Fig. 8. By 
contrast, stars in Regions D and E have much larger flux 
density excesses at 60 /zm, which make them detectable 
at much lower total flux density levels. Stars in Regions B 
and G also differ markedly from those in Regions D and 
E as far as the [60] — [100] index is concerned (Fig. 9; see 
also Fig. 6): the resolved envelopes in Regions D and E go 
along with large 100 /im excesses, suggestive of cool dust 
in detached envelopes, contrary to the situation prevailing 
in Regions B and C (the only exception being T Sgr, which 
appears as a border case between C and D) . In the case of 
RZ Sgr, which is extended in both the 60 /um (YPK) and 
100 um bands (according to the IRAS Small Scale Struc- 
ture Catalogue 1985), an optical nebula has even been re- 
ported by Whitelock (1994). 

In Region E, Ori and S929 do not follow the general 
trend. In those cases, there may be a confusing backgroimd 
source responsible for the strongly asymmetric 60 /zm pro- 
file (see Fig. 10), as the 100 /xm profile is offset by about 
1' in the direction of the 60 /im asymmetry. The extension 
observed for S929 may be real; its [60] — [100] color index 
changed by only 2% between the two IRAS observations, 
while the 60 /im flux density changed by 20%. This color 
stability is observed in other variable IRAS sources like x 
Cyg (see Appendix A), but would probably not be pre- 
served if the excess flux density were due to a background 
source. 

Finally, the flux density excess in Regions D and E 
is associated with a widening of the whole 60 /im (and 
sometimes 100 /xm) source profile, whereas the excess for 
stars in Regions B and C is caused by a weak plateau at the 
base of the profile, as apparent on Fig. 10. For the resolved 
shells in Regions D and E, the full widths at the 25% and 
50% flux density levels (Table 2) are indeed much larger 
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Fig. 8. Flux ratio F^/Fp vs F(2.2). Symbols are as in Fig. 7, 
except that large squares now denote stars fulfilling simultane- 
ously our two criteria for extended envelopes. Stars have been 
separated according to Regions B, C (upper panel) and D, E 
(lower panel) 



than those expected for a point source. HD 191630 and RZ 
Sgr are listed in the IRAS Small Scale Structure Catalogue 
as extended at 60 /um and 100 /um, respectively. In fact, 
this distinctive feature of the resolved shells in Regions 
D and E has been used to include in Table 2 two stars 
(BI And and AA Cam) with wide profiles, despite quality 
factors QF < 5 that would normally not qualify them. 
However, these are distant stars with small 60 jjm flux 
densities, so that the F^/Fp ratio cannot be determined 
accurately (and is therefore not listed in Table 2). 

Stars with resolved shells are represented as black dots 
in the {K-[12], [25] - [60]) diagram (Fig. 11). This figure 
illustrates the relative contribution of the extended shell 
to the [25] — [60] index: a line segment joins the [25] — 
[60] indices computed from the 'zero- crossing' flux den- 
sity Fz at 60 fim (measuring the combined contributions 
of the star and its resolved shell) and from the 'template' 
flux density Ft, which is assumed to be a rough measure 
of the stellar contribution alone (represented by a black 
dot in Fig. 11). Although that assumption is certainly a 
very rough one, it is not totally unreasonable as stars from 
Region E now move down, and most reach Regions B, C 
and D when adopting F^ instead of F^ to represent the 
60 /xm photospheric flitx density. Some stars, however, do 
not quite leave Region E, presumably because the point 
source fitting yielding Ft does not entirely remove the con- 
tribution of the resolved shell in those cases. 
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Fig. 9. Same as Fig. 8 for F^/Fp vs [60] - [100] 
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Fig. 10. The 60 ^m co-added, spline-interpolated scans 
(dashed lines) for all stars with a possibly resolved shell in 
Regions D - E, and for selected cases in Regions B - C (see 
Table 2) . The dotted line is the template 60 /im profile as pro- 
vided by IPAC. Note how larger the deviations from the tem- 
plate are in Regions D - E as compared to Regions B - C. In 
the latter case, the deviation of the observed profile from the 
template profile is due to a weak extended tail. Such detections 
are therefore vulnerable to detector hysteresis (see text) 
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Fig. 11. Position of the sources with resolved shells (filled 

circles) in the the {K - [12], [25] - [60]) diagram. The filled 
circles represent the color computed from the 'template' flux 
density Ft at 60 /um, i.e. they roughly correspond to the flux 
density of the star alone. The upper end of the vertical segment 
is located at the color computed from the 'zero-crossing' flux 
density Fz at 60 /um, i.e. the shell + star color 
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Fig. 12. Location of the maser sources in the {K — [12], [25] — 
[60]) diagram. Positive detections correspond to filled symbols, 
non-detections to open symbols, and stars not probed for maser 
emission to crosses. Squares stand for SiO masers, triangles for 
H2O masers and circles for OH masers 



Especially interesting is the fact that Y Lyn and OP 
Her move from Regions C and D, respectively, to Region 
B, which is well in line with their small ZrO index (see 
Fig. 3 where they appear as outliers). With this adjust- 
ment, all three SRc variables in our sample (RS Cnc, Y 
Lyn and T Get) now belong to Region B, and moreover 
have resolved shells (tentative in the case of T Cot). As 
noted by Young et al. (1993a,b) and Habing (1996), re- 
solved shells are a common property of semi-regular vari- 
able stars. 

Finally, it should be mentioned that several of the stars 
with a shell resolved by IRAS turn out to have an extended 
CO shell as well, as derived in Sect. 5.3 from the modelling 
of the CO data. The inferred radius of the CO shell (Ta- 
ble 6) is larger than 10" for x Cyg (18"), W Aql (23"), 

Gru (15") and FU Mon (60"). 

3. Masers in S stcirs 

This section presents data collected from the literature 
on SiO, OH or H2O maser emission from S stars (relying 
on the compilation of Benson et al. 1990 for the earlier 
literature). These data are summarized in Table 3. The 
detection of SiO, OH or H2O maser emission is a clear 
indication that the circumstellar shell is oxygen-rich (see 
e.g. Deguchi et al. 1989, Lewis 1996), thus providing an 



important constraint on the chemical nature of the dust 
grains (see Sect. 4). In Table 3, 'Y' or 'N' in the SiO, OH 
or H2O columns means that the corresponding maser is or 
is not present, respectively, whereas a dash indicates that 
the maser has not been searched for in a given star. 

The masers which are found in S stars are mainly SiO 

masers forming in the densest part of the circumstellar 
envelope near the stellar photosphere. OH maser emis- 
sion has been detected in RS Cnc by Rudnitskij (1976) 
but not by Kolena & Pataki (1977), though the latter au- 
thors probed a different transition from that (1667 MHz) 
detected by Rudnitskij (1976). 

The S stars with SiO maser emission lie in the lower 
part of Region C, with some overlap with Regions B and D 
(Fig. 12). As expected, the region occupied by the masers 
exactly matches the region delineated by the stars showing 
the 9.7 /Ltm silicate feature in emission (IRAS LRS class 
E; Fig. 5). Indeed, 6 of the 8 SiO masers with an available 
LRS class are of class E, the only exceptions being R Lyn 
and EP Vul (class F). On the other hand, no maser emis- 
sion is observed for stars in Regions D and E (with the 
exception of EP Vul, but that star lies very close to the 
boundary with Region C). 
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4. A simple model of dusty circumstellar shells 

4.I. Ingredients of the model 

In order to relate the diversity of IR color indices observed 
among S stars to the underlying physical parameters, syn- 
thetic IR color indices of a star embedded in a circumstel- 
lar shell have been computed for various input parameters. 
In our simple model, the star is assumed to radiate as a 
black body at a temperature Toff- The mass- losing star 
is surrounded by a spherically-symmetric dust shell ex- 
tending from Tin to Tout) with density decreasing as r~^, 
r being the distance from the central star. This is equiv- 
alent to assuming a steady mass-loss rate M at constant 
outflow speed. The inner radius of the dust shell must be 
larger than or equal to the radius where grains start con- 
densing (i.e., to the radius where the shell temperature 
drops below 1300 K for silicates, or 1500 K for graphite 
and amorphous carbon). The shell outer radius is chosen 
such that Tout /''in = 10^, which ensures that the color 
indices of the shell have reached an asymptotic value. A 
roughly logarithmic radial mesh is defined in the dust en- 
velope so that each shell is optically thin. In each shell, 
the grains are assumed to be in thermal equilibrium, so 
that the energy absorbed by the grains exactly balances 
the energy re-emitted. At the inner boundary, the radi- 
ation field is that of a black body of temperature Teff. 
The model IRAS flux densities are calculated from the 
emergent spectrum by convolving it with the IRAS fil- 
ter bandpass (IRAS Explanatory Supplement, 1988). The 
main shortcoming of the code is the neglect of the scat- 
tering contribution, since only absorption is taken into 
account. 

Three types of dust grains have been considered: sili- 
cates (with a specific mass of 3.5 g cm~^), graphite (with 
a specific mass of 2.25 g cm^-^) and amorphous carbon 
(with a specific mass of 1.85 g cm~^). The grain radius is 
0.2 in all cases. The absorption coeflBcients as a func- 
tion of wavelength are taken from Draine & Lee (1984) and 
Draine (1985) for silicates and graphite. For amorphous 
carbon, the absorption coefficients have been generated 
with the usual Mie formulae using the optical constants 
provided by Rouleau & Martin (1991). At wavelengths 
A > 50 /Lim, the spectral index of the emissivity coeflacient 
has been taken equal to —2 for graphite grains, and to 
—1.5 for amorphous-carbon and silicate grains (Ivezic & 
Elitzur 1995). 

4-2. Synthetic color indices 

The results of this model are presented in Fig. 13 for dust 

shells made of either silicate grains, graphite grains or 
amorphous carbon grains. Constant dust mass loss rates 
of 10-12^ 10-1°, 10-^ 10-s and 10"^ M© y-^ with a wind 
velocity of 14 km have been adopted. The shell in- 
ner radius is set by the grain-condensation temperature, 
so that the dust shells computed in Fig. 13 are not 'de- 



tached' (in the sense of Willems & de Jong 1988). The 
central star has been assigned effective temperatures Tes 
of 4000 K (sohd line) and 3000 K (dashed line), and a lu- 
minosity L = 5000 Lq ; these parameters turn out to have 
little impact on the shell colors. 

In the ([12] - [25], [25] - [60]) and (K - [12], [25] 

— [60]) diagrams, the silicate track is quite distinct from 
the graphite or amorphous-carbon tracks. The graphite 
and amorphous-carbon tracks go directly from Region I 
to upper VII (i.e., from A to D and upper C), whereas 
the silicate track goes from I to II and Ilia (i.e., from A 
to B and lower C). These differences observed in Fig. 13 
between carbon- and oxygen-rich shells must be related 
to the different emissivities of silicate and carbonaceous 
grains in the IRAS bands, as discussed by Ivezic & Elitzur 
(1995). 

The tracks for carbonaceous and silicate grains pre- 
dicted by this simple model outline the segregation ob- 
served in the color-color diagrams between S stars with 
an oxygen-rich shell, as indicated by the 9.7 /zm silicate 
feature (IRAS LRS class E), and S stars with featureless 
IR spectra (IRAS LRS class S 'stellar' or F 'feature- 
less'). It is important to note here that featureless spec- 
tra are indeed predicted for carbonaceous grains (Ivezic & 
Elitzur 1995). 

The data on mascr emission collected in Sect. 3 are 
compatible with this segregation. Those S stars which are 
SiO or OH masers (Table 3) must have an oxygen-rich 
circumstellar environment, and lie indeed close to the sili- 
cate track (Fig. 12). On the contrary, no SiO, OH or H2O 
maser emission has been detected for S stars in Regions D 
and upper C, along the tracks corresponding to carbona- 
ceous grains. Those S stars might therefore possibly be 
surrounded by C-rich circumstellar shells, especially since 
they occupy the same region of the color-color diagram as 
the optical carbon stars (Chan & Kwok 1988). 

However, the above picture is not totally satisfactory, 
since (i) the S stars with silicate emission are actually 
located in between the silicate and carbonaceous tracks, 
(ii) many stars in Regions D and E have large 60 /xm ex- 
cesses that cannot be reproduced by the tracks displayed 
in Fig. 13, (iii) at least one star (R Gem) moves from the 
region of silicate dust (lower C) into the region of car- 
bonaceous dust (upper C) during its variability cycle, and 
(iv) the photosphere of S stars is oxygen-rich, so that their 
circumstellar shell may be expected to be oxygen-rich as 
well. The first two mismatches are in fact not specific to 
S stars, and possible solutions were already suggested by 
Ivezic & Elitzur (1995). They include either (a) invoking a 
mixture of silicate and carbonaceous grains (see however 
the discussion of Sect. 3 on maser emission), (b) decreas- 
ing the spectral index of the emissivity at long wavelengths 
to values smaller than —1.5, or (c) considering detached 
envelopes in the sense advocated by Willems & de Jong 
(1988). 
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Fig. 13. Left panel: The ([12] — [25], [25] — [60]) colors predicted for circumstellar shells (with rout/nn ~ 10*) made of silicate 
grains, graphite grains or amorphous-carbon grains, surrounding a star with L = 5000 Lq and Tcs— 4000 K (solid lines) or 
3000 K (dashed lines). The diamonds along the curves correspond to dust mass loss rates of 10"^^ 10"^°, 10"^ 10"* and 
10~^ M0 (from left to right), with a wind velocity of 14 km s~^. Black bodies fall along the central straight line labelled 
BB. Right panel: same as left for {K - [12], [25] - [60]) 



As discussed in Sect. 2.4, the S stars located in Re- 
gions D and E have very specific properties that may help 
to identify the origin of their large 60 /im excess. First, 
the prototypical SC stars GP Ori and FU Mon located in 
Region E have a C/0 ratio equal to unity to within 1% 
(Dominy et al. 1986). Similarly, TT Cen is a rare CS star 
which exhibits ZrO bands at some times and C2 bands at 
other times (Stephenson 1973). According to Stephenson 
(1973), these variations are probably caused by temper- 
ature changes in an atmosphere with a C/0 ratio very 
close to unity, or perhaps even to a secular change in the 
atmospheric C/0 ratio. Most of the stars in Regions D 
and E arc actually SC stars or at least S stars with a large 
C/0 spectral index. Finally, several stars in Region E have 
very extended shells that are resolved at 60 yum and, in one 
case (RZ Sgr), visible in the optical (Whitelock 1994). All 
these properties point towards these stars being in a very 
rare and short-lived evolutionary phase. Based on this ev- 
idence, we suggest that the large 60 fim excess of SC stars 



populating Regions D and E finds a natural explanation 
in the much debated concept of interrupted mass loss first 
proposed by Willems & de Jong (1988). The cessation of 
mass loss when C/0 gets close to unity (because all C 
and O atoms are then locked into the CO molecule in- 
stead of being involved in dust-forming molecules) causes 
the dust shell to detach from its parent star, and to cool 
down as it expands into the interstellar medium without 
being replenished at its inner side. As shown by Willems 
& de Jong (1988) and Chan & Kwok (1988), the colors 
of the dust shell then describe a counter-clockwise loop in 
the color-color diagram, starting from the region of stars 
with silicate emission and ending close to the photospheric 
point (Region A) when the shell has dissolved into the in- 
terstellar medium, after passing through Region E. The 
mass loss resumes when C/0 reaches values above unity, 
and the star then enters the region of heavily-obscured 
infrared carbon stars (Chan & Kwok 1988; lower VII in 
Fig. 13). The SC stars found in that region (RZ Peg, UY 
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Cen and BH Cru, a sister case of TT Cen uncovered by 
Lloyd Evans 1985) may actually be on the lower part of 
that loop. The application of the Willcnis & dc Jong sce- 
nario (implying that the loop described in the color-color 
diagram corresponds to a brief evolutionary phase) to the 
rare, supposedly short-lived SC phase would thus not face 
the difficulty of inconsistent time scales generally used to 
argue against it (e.g., Zuckerman & Maddalena 1989). Fur- 
ther support for this idea comes from the peculiar CO 
line profiles observed for many of the stars populating Re- 
gion E (FU Mon: see Fig. 15 below; DK Vul, RZ Sgr, TT 
Cen and UY Cen: Sahai & Liechti 1995), as discussed in 
Sect. 6.3. 

5. Molecular line data and mass loss rates 

5.1. New CO data 

Observations of the C0(2— 1) line at 230 GHz and of the 
C0(3— 2) line at 345 GHz were made in January 1996 and 
January 1997 with the 10.4 m telescope of the Caltcch 
Submillimeter Observatory (CSO) on Mauna Kea, Hawaii. 
The CSO is equipped with SIS junction receivers cooled to 
liquid helium temperatures. The effective single-sideband 
system temperatures for these observations, including the 
effects of atmospheric emission and absorption, were about 
500 K and 800 K at 230 and 345 GHz, the telescope half- 
power beamwidths were respectively 30" and 20" and the 
main-beam efficiencies 76% and 65%. 

The spectra were observed using two 1024 channel 
acousto-optic spectrographs (AOS) simultaneously. The 
first has a total bandwidth of ~ 500 MHz (~ 800 km s"^ 
at 230 GHz) and a velocity resolution of ^ 1 km s""'^. The 
second has a bandwidth of ~50 MHz and a velocity resolu- 
tion of 0.1 km s~^. The observations were made by chop- 
ping between the star position and an adjacent sky posi- 
tion, offset 90" in azimuth, at a rate of 1 Hz, and consisted 
of pairs of chopped observations with the source placed 
alternately in each beam. The spectral baselines resulting 
from this procedure are linear to within the r.m.s. noise. 
The telescope pointing errors were measured by mapping 
the spectral line emission from a nearby CO bright star 
before each observation was made, and the pointing accu- 
racy is better than ~ 5" for all of the observations. 

The temperature scale and atmospheric opacity were 
measured by chopping against a hot (room temperature) 
load. The line temperature was corrected for the main- 
beam efficiency, and the resulting scale is the Raylcigh- 
Jeans equivalent main beam brightness temperature Tmb , 
i.e. that measured by a perfect 10.4 m antenna above the 
atmosphere. The spectrometer frequency was calibrated 
using an internally generated frequency comb, and the ve- 
locity scale is corrected to the Local Standard of Rest 
(LSR). 

Twelve S stars were observed at the CSO (Table 4) , 
nine in the C0(2-l) line and three in the CO(3-2) line, 
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Fig. 14. C0(2-l) line profiles of four S stars observed at CSO. 
The abscissa is velocity with respect to the LSR and the or- 
dinate main beam brightness temperature. The profiles for 
GCGSS 212 and 589 are observed with a velocity resolution 
of ~ 0.1 km s"\ while GCGSS 283 and GCGSS 704 are ob- 
served with 1 km s~^ resolution. The horizontal bars show the 
velocity range observed in optical spectra (Table 4). There is 
no optical velocity available for GCGSS 704 

and emission was detected from four (Table 5). The line 

profiles for the detected stars are shown in Fig. 14. 

Table 4 lists the stars that were observed: the GCGSS 
number and the variable star name are in columns 1 and 
2, respectively. Next is the observed position: we used po- 
sitions accurate to ~ 1" from the HST Guide Star Catalog 
and other sources (see Chen et al. 1995). Columns 5-9 list 
the galactic longitude and latitude, the spectral type from 
the GCGSS, the variable type and the period from the 
GCVS. Column 10 gives the stellar radial velocity with 
respect to the LSR measured from optical spectra; since 
the radial velocities of red giants vary as they pulsate. Ta- 
ble 4 lists the range of reported radial velocities. Finally, 
column 1 1 gives the channel-to-channel r.m.s. noise in the 
500 MHz AOS. 

Table 5 gives the line parameters of the four detected 
stars: the C0(2 1) line flux in K x km s~^, the peak line 
temperature, the central velocity and the half-width of 
the line at zero power, Vf,, which gives the terminal wind 
outflow speed. These quantities are determined by fitting 
a parabolic line model to the data. The agreement between 
the optical and CO radial velocities is good. 

Seven of the stars in Table 4, GCGSS 89, 117 (CP Ori), 
422 (NQ Pup) 626 (FM Hya), 704 (Z Ant), 796 (HR 4755) 
and 803 (S UMa) have not previously been observed. We 
detect one of these, Z Ant. GCGSS 816 (UY Cen) is weakly 
detected by Sahai & Liechti (1995, SL95) with good agree- 
ment between the CO and optical radial velocities. This 
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Fig. 15. Broad-band C0(2-l) spectrum in the direction of 
GCGSS 212 (FU Mon). The horizontal bar shows the range of 
the observed optical velocities 



star was not detected in the present observations, but our 
sensitivity is lower. The detection of circumstellar CO for 
this star is of particular interest, given its rare SC spectral 
type (sec Sect. 6.3). GCGSS 149 (NO Aur) and 796 (HR 
4755) were previously observed by SL95 and by Bieging & 
Latter (1994, BL94); like the present observations, these 
did not detect CO emission. 

RS Cnc has been observed by many authors (see Loup 
et al. 1993, for example). Margulis et al. (1990) point out 
that the CO line profile for this star, as for several oth- 
ers, more closely resembles a triangle or gaussian in shape 
than the parabolic profile typical of circumstellar winds. 
The high velocity resolution observations in Fig. 14 show 
that the line profile actually consists of two parabolic com- 
ponents of different widths centered at the same velocity. 
The parameters for these profiles, estimated by eye, are 
given in Table 5. This line shape may indicate the pres- 
ence of two molecular winds. Such line profiles have been 
seen for several other stars (e.g. Margulis et al. 1990, SL95, 
Kahane & Jura 1996, Knapp et al. 1997a) and may be 
quite common. 

GCGSS 283, R Lyn, has previously been detected by 
BL94, and the data in Table 5 are in good agreement. FU 
Mon (GCGSS 212) was detected by SL95, who observe 
two narrow features at —44 and -1-14 km s^^ which they 
attribute to the blue- and red-shifted components of an 
expanding circumstellar shell. However, the optical radial 
velocity of the star is about —43 km (Table 4), and 
further, the narrow component at -|-14 km is likely 
to be interstellar, as shown by our observations with the 
500 MHz AOS in Fig. 15. We conclude that the emission 
at —44 km s~^ (best seen on the high resolution profile 
in Fig. 14) is from the circumstellar envelope. The out- 



flow speed, 3 km s ^, is very low, but similar low values 
are found for some other Mira and semi-regular variables 
(Young 1995; Kerschbaum et al. 1996). 

5. 2. CO data from the literature 

Table 6 summarizes (under the header 'Observations') CO 
millimeter wavelength observations of S stars published 
since 1990, including those in the present paper. Data ob- 
served prior to 1990 can be found in Loup et al. (1993). 
The stars in Table 6 are grouped according to their loca- 
tion in the {K — [12], [25] — [60]) color color diagram (see 
Sect. 2.3) and arranged in order of right ascension within 
these groups. Several observations are not listed in Ta- 
ble 6 because they were made at positions which are too 
discrepant from the optical position; the stars are GCGSS 
133, GZ Peg, and T Cam. 

Table 6 gives the IRAS name (an asterisk before the 
IRAS name refers to a note at the end of the table), the 
star name and the results of CO observations of the star: 
the line observed; the telescope half power beamwidth in 
arcseconds; the channel-to-channel r.m.s. noise in K; the 
integrated line brightness in K x km s~^; the peak bright- 
ness temperature in K; the central velocity Vc with respect 
to the LSR; the wind outflow speed 14; and the reference. 
All temperatures are expressed in main-beam brightness 
temperature. Dashes for any of these quantities mean ei- 
ther that no emission was detected from the star or that 
the quantity in question was not quoted in the paper. Ta- 
ble 6 contains observations of 56 stars, with 35 detections. 

To first order, the peak brightness temperatures and 
the integrated CO line intensities should scale inversely 
with the square of the telescope beamwidth, since the CO 
lines are usually fairly optically thick and the envelopes in 
general smaller than the beam. As Table 6 shows, this is 
roughly the case, and the agreement among the observa- 
tions is in general good, with no serious discrepancies. 

Figure 16 shows the histogram of the stellar systemic 
radial velocities (with respect to the LSR). The values 
adopted for stars with multiple observations are straight 
averages of the individual values. The mean velocity for S 
stars detected in CO is < V > = —9.1 ± 5.0 km and 
the radial- velocity dispersion is a = 25.9 ± 4.5 km s~^. 
This value refers to a sample of intrinsic S stars, since no 
extrinsic S stars have been detected in CO. This dispersion 
is typical of a young-disk population (Mihalas & Binney 
1981). Such a population should have a scale- height above 
the galactic plane of about 200 pc, in agreement with the 
results of Van Eck et al. (1997). 

Figure 17 shows the histogram of the wind termi- 
nal velocity Vc compared with the distributions for three 
other sets of molecular line observations; those for nearby 
oxygen-rich Mira variables (Young 1995), for semi-regular 
(SRa and SRb) variables (Kerschbaum et al. 1996) and for 
carbon stars (Olofsson et al. 1993). The outfiow speeds 14 
for the S stars are taken from Table 6. We used average 
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Fig. 16. Distribution of radial velocities (with respect to the 
LSR) from CO observations of intrinsic S stars 



values for stars with multiple observations except for RZ 
Sgr (20120-4433), for which CO(l-O) and (2-1) obser- 
vations give discrepant values (14 and 8.8 km s^^; SL95). 
We use the velocity derived from the CO (2-1) observation 
since this has a much higher signal-to-noise ratio, but note 
that the CO(l-O) line may really be broader; the larger 
telescope beam at this wavelength could be detecting gas 
at larger distances from the extended envelope of this star, 
which could have a larger outflow velocity. 

Figure 17 shows that oxygen-rich Miras have the small- 
est outflow velocities (median 6.5 kms~^, largest value 
12.7 km s~^), while those of the SRVs cover a similar 
range (median 8.0 kms""*^, largest value 15.6 km s^^). 
Carbon stars have the largest outflow velocities (median 
12.0 km s~^, largest value 33.2 km s~^) while as expected 
S stars are intermediate (median 8.5 km s~^, largest value 
24.7 km s~^). The largest Ve in our sample of S stars is 
observed for the CS star TT Cen, a border case between S 
and C stars (see Sect. 2.1). Figure 18 shows that, among 
Mira S stars, the outflow velocity correlates well with the 
period of the photometric variations, a result already dis- 
cussed by Heske (1990) and by Olofsson et al. (1993). Jura 
(1988) finds an almost identical dependence on period for 
the ratio of flux densities at 25 and 2.2 /xm. 



5.3. Mass loss rates 

The CO data from Table 6 were used to calculate mass 
loss rates for the detected S stars and upper limits for 
the non-detected stars. The CO lines were modeled us- 
ing a code based on that of Morris (1980) which assumes 
spherically-symmetric mass loss at a constant rate and 
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Fig. 17. Normalized histogram of terminal wind outflow 

speeds of (a) oxygen-rich Miras (Young 1995) and somircgu- 
lar variables (Kerschbaum et al. 1996) and (b) S stars (present 
work) and carbon stars (Olofsson et al. 1993) 



constant outflow speed^ with excitation by collisions and 
by infrared photons at 4.6 /im. The envelope outer radius 
was taken to be that at which the CO is photodestroyed 
by the diffuse interstellar radiation field using the calcu- 
lations of Mamon et al. (1988). The details are given by 
Knapp et al. (1997b). The relative abundance of CO to 
H2 was assumed to be 6.5 x 10"'* for all stars (Lambert et 
al. 1986; Smith & Lambert 1990). 

The infrared radiation field was approximated as that 
of a black body of temperature 2500 K and radius 2.5 10"'^^ 
cm. Models show that the CO line strength is only weakly 
dependent on the radiation field, so this simplifying as- 
sumption is unlikely to produce an uncertainty of more 
than 20% in the derived mass loss rates. 

The distances are derived by adopting absolute mag- 
nitudes which depend on location in the IR color-color 
diagram. For stars in Region A (photospheric colors) we 
assume Mboi = —2, i.e. one magnitude below the RGB 
tip for solar-metallicity stars with M = 1 Mq (Schaller et 
al. 1992). The bolometric correction Mboi — Mk is derived 
from the apparent bolometric fiux for BD Cam (HR 1105) 
obtained by integrating the flux densities corresponding 
to the UBVRIJKL magnitudes from Lee (1970) plus the 
IRAS flux densities. The corresponding absolute K mag- 



Thc mass loss rates estimated in the present paper might 
therefore be inadequate for those stars in Regions D and E 
with a resolved - possibly detached - shell as observed at 60 
/um 
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Fig. 18. Terminal wind outflow speed of S stars measured from 

CO millimeter wavelength emission lines versus period of pho- 
tometric variations (from the GCVS). Circles: Mira variables 
(open circles are for stars for which there is only one measure- 
ment of Vc); crosses: semi- regular variables. The points corre- 
sponding to the two outflow speeds for RS Cnc are connected 
by a vertical line. The data for DY Gem, 06331+1414 (SRa; 
P = 1145 d), are not included in the figure 



type SC, the choice Mk = — 8.1 appears to be a reasonable 
one. 

Wc modeled the wind from RS Cnc as two separate 
components, fit only to the C0(2-l) observations de- 
scribed in the previous section. 

The upper limits to the mass loss rate of S stars not 
detected in CO were calculated using the median outflow 
speed of 8.5 km s"^ found for the detected sample and 
assuming that we can detect a line of brightness tempera- 
ture three times the r.m.s. noise (examination of the data 
in Table 6 suggests that this is reasonable). 

The best fit mass loss rates were found by calculating 
the model peak antenna temperature and integrated line 
intensities for a given input mass loss rate and comparing 
it with the observations. The mass loss rates were adjusted 
until reasonable agreement with all of the observations 
was found, and were corrected for the mass of helium (see 
Knapp & Morris 1985). 

The results are listed in the rightmost columns of Ta- 
ble 6 (under the header 'Model'), which list the distance, 
the mean outflow speed (the same for all lines; the value of 
8.5 km s~^ used for calculating the upper limits is given 
in parentheses), the CO photodissociation radius in cm 
and in arcseconds, and the predicted CO peak line tem- 
perature and integrated CO line intensity for each of the 
observed lines. Comparison between the calculations and 
observations shows that uncertainties in the observations 
introduce about a factor of two uncertainty into the mass 
loss rate. 



nitude is —4.6. Other authors, e.g. Jura (1988) have used 
Mk = —8.1 for all S stars; this absolute magnitude is 

derived from carbon stars in the solar neighborhood and 
AGB stars in the Magellanic Clouds. 

For stars in the other regions, which are supposedly 
more evolved than are those in Region A, we adopt 
Mk = —8.1. That value of Mk yields distance moduli 
that are consistent with direct determinations when avail- 
able. For the composite system tt^ Gru (S5,7e + GOV), 
Ake & Johnson (1992) derive a distance modulus of 6.0 
from a fit to the UV spectrum, corresponding to a dis- 
tance of 160 pc, identical to the value derived from the 
K magnitude (Table 1). For T Sgr, a distance of 1000 pc 
(as compared to 810 pc from the K distance modulus) is 
derived by Culver & lanna (1975) from the spectral type 
F3IV assigned to its companion. For x Cyg, a distance 
of 136 pc is obtained from the distance modulus in the 
K band, consistent with that (106 pc) derived from the 
Hipparcos parallax (tt = 9.43 ± 1.36 mas; van Leeuwen et 
al. 1997). 

The evolutionary status of stars in Regions D and E is 
imclear, as is their relationship to stars in other regions of 
the color-color diagram, and so their absolute magnitudes 
are uncertain. Since these are generally stars of spectral 



6. Discussion 

6.1. Comparison with mass-loss rates of Oxygen and Car- 
bon stars 

Several authors (e.g. Netzer & Elitzur 1993; Young 1995) 
have shown that AGB stars with higher mass loss rates 
have higher outflow velocities. The corresponding relation- 
ship for S stars is shown in Fig. 19 and compared with 
the results for oxygen-rich Miras (Young 1995) and car- 
bon stars (Olofsson et al. 1993). The relationship between 
Vc and M found for oxygen stars by Young (1995) is also 
shown. 

It is not clear at present whether the relationship be- 
tween Vc and M demonstrated by these data is real, in the 
sense that it has a physical origin, or is due to selection ef- 
fects. For the present purpose, we will simply treat it as a 
convenient way to display the data and compare samples. 
The mass loss rates in Fig. 19 have been calculated as- 
suming a CO abundance of / = n(C0)/n(H2) = 3 x 10~^ 
for oxygen stars (Yoimg 1995), 6.5 x lO^** for S stars and 
9 X 10"'' for carbon stars (Olofsson et al. 1993). The S and 
oxygen stars show essentially the same correspondence be- 
tween M and Vc (though there is a lot more scatter in the 
S star data), giving some confidence in the assumed value 
of / and the resulting mass loss rates. The apparently 
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Fig. 19. Wind outflow speed T4 versus mass loss rate for: open 
circles: oxygen-rich Miras (Young 1995); filled circles: S stars 
(current work); crosses: carbon stars (Olofsson et al. 1993). 
The relationship for oxygen Miras found by Young (1995) is 
also shown 

higher values of Vq for a given mass loss rate for carbon 
stars (see also Fig. 17) will be discussed elsewhere. 

6.2. Comparison of mass loss rates and color excesses 

The color excess can also be used to give the mass loss rate 
if the gas-to-dust ratio and the outflow speed are known. 
This is illustrated in Fig. 20, which shows the envelope 
density M /Vg versus the 12/;im/2.2/xm color. At low mass 
loss rates (M < 10~^ M0 y~^) the broad-band colors 
arc dominated by the colors of the photosphere, as ex- 
pected, while the stars with higher mass loss rates show 
a strong correspondence between the mass loss rate and 
the 12 /im color, showing that the gas to dust ratio in 
these envelopes is roughly constant. The relationship in 
Fig. 20 is in agreement with that for oxygen stars shown 
by Habing (1996). 

The two stars TT Cen and RZ Sgr in Region E 
markedly depart from this relationship, however. These 
are probably in a transitory phase of evolution; TT Cen is 
a rare CS star where ZrO bands seem to have disappeared 
while C2 bands appeared (see Sect. 2.1 and Stephenson 
1973), while RZ Sgr is surrounded by an optical (White- 
lock 1994) and IR (YPK) nebula. 

Stars with roughly photospheric 12/im/2.2/im colors 
show a wide range in envelope densities. This large scatter 
may be due to the imperfect coupling between dust and 
gas at these low densities, which sets a lower limit to the 



Fig. 20. Envelope density M /Vc (measured in M0 y~^ and 
km s^^) versus the ratio of 12/im to 2.2/im flux densities. The 
vertical dotted line shows the photospheric color for a star of 
temperature 3000 K. Open symbols: stars in Regions A, B and 
C of the IR color-color diagram (Fig. 3); flUed squares: stars 
in Region D; filled circles: stars in Region E. The inverted 
triangles show upper limits 

mass loss rate for a radiation-pressure driven wind (cf. 
Netzer & Elitzur 1993; SL95). Empirically, this hmit is a 
few 10-^ Mq y-i (Fig. 20 and Table 6). 

6.3. Mass loss rates and the IRAS color-color diagram 

Figure 21 presents the mass loss rates of S stars as a func- 
tion of their location in the (i^T — [12], [25] — [60]) diagram. 
None of the stars in Region A (stars with photospheric col- 
ors) has detectable circumstellar CO emission, with limits 
on the mass loss rates of < 6 x 10~^ Mq and envelope 
densities well below those of the detected stars (Fig. 20). 
This confirms the lower limit at which a star can lose mass 
by a radiation-pressure driven wind estimated by Netzer 
& Elitzur (1993). The low mass- loss rates inferred for ex- 
trinsic S stars in Region A also confirm an earlier sugges- 
tion (Paper I) that these stars are less evolved than the 
intrinsic S stars populating the other Regions of the IR 
color-color diagram. 

Stars in Regions B and D have moderate mass loss 
rates, in the range 5 10^*^ to 5 10^^ M0 y^^. One star 
from Region B (RS Cnc) has been found to exhibit a dou- 
ble wind (Table 6). Most of the observed stars, and most of 
the detections, lie in Region C, which contains stars with 
moderately optically thick circumstellar envelopes, likely 
containing silicate dust. CO emission is detected from 20 
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Fig. 21. Mass loss rates of S stars (as measured from CO) 
in the {K - [12], [25] - [60]) diagram. The diameter of the 
circle is proportional to the mass loss rate, as labeled. Inverted 
triangles correspond to upper limits 



Fig. 22. The outflow velocity as measured from CO for S stars 

in the {K — [12], [25] — [60]) diagram. The diameter of the circle 
is proportional to the outflow velocity. Compare to Fig. 10 of 
Olofsson et al. (1993) for C stars 



of these 23 stars, and the mass loss rates are typically 
larger than several xlO^^ M0 y^^. The undetected stars 
have upper limits greater than this value, so that the data 
are consistent with the conclusion that all stars in this re- 
gion lose mass at a rate larger than several 10^^ Mq y^^. 

The stars in Region D, which show roughly photo- 
spheric 12/um/2.2/um colors despite large [25] — [60] in- 
dices, generally have low mass- loss rates (like Region B). 
Stars in Region E have mixed properties. Some, like TT 
Cen and RZ Sgr, lose mass at a very large rate (several 
10"^ Mq y-1), while others like FU Mon lose mass at 
a more moderate rate (a few 10"'' Mq y~^). Figure 22 
presents the variation of the wind velocity across the 
(K — [12], [25] — [60]) diagram, and shows that stars in 
Region E also have a wide range of wind velocities, from 
very low (FU Mon: 2.8 km s~^) to very large (TT Cen: 
24.7 kms-i). 

FU Mon also has a resolved IR envelope (Sect. 2.4), 
and this fact taken together with a low mass loss rate and 
a small outflow velocity suggest that mass loss has just re- 
sumed in that star. Olofsson et al. (1990; 1993) detected 
several carbon stars (S Set, U Ant, TT Cyg) with a de- 
tached shell and a double wind that would also fall in our 
Region E. The older, detached shell was produced by a 
massive, fast wind, whereas the recent shell is caused by a 
slow (5 km s""'^), low mass loss rate wind. Since FU Mon 
is an SC star with a C/0 ratio close to unity (Dominy et 
al. 1986), the mass loss may have come to a halt when the 



the C/0 ratio approached unity, as already suggested by 
Willems & de Jong (1988) and Chan & Kwok (1988). This 
scenario may well hold for all SC stars, even though some, 
like UY Cen, are in fact located in Region C. As argued 
in Sect. 4.2, that star may be at the end of the loop in 
the IR color-color diagram. The horn-shaped CO line pro- 
file of UY Cen observed by SL95, indicative of a detached 
shell and fossil mass loss, supports this idea. Finally, other 
stars lying in Region E like TT Cen, RZ Sgr and DK Vul 
have peculiar CO line profiles (narrow central feature su- 
perimposed on a broader less intense feature) suggesting 
that they have multiple winds. 



6.4- Mass loss rates and binaries 

Several intrinsic S stars in the sample considered in this 
paper are binaries with main sequence companions, as re- 
vealed by the composite nature of their spectrum at min- 
imum light (WY Cas, W Aql, T Sgr; Hcrbig 1965; Culver 
& lanna 1975). A composite spectrum is also suspected for 
S Lyr from its shallow lightcurve (Merrill 1956). The star 
TT^ Gru has a close GOV visual companion (Feast 1953). W 
Aql, WY Cas and S Lyr are surely among the S stars with 
the largest mass loss rates, while T Sgr and tt^ Gru have 
mass loss rates close to the average for stars in Region C. 
Clearly, a definite conclusion as to whether binarity can 
indeed reinforce the mass loss rate, e.g. by the companion- 
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reinforced attrition process (Eggleton 1986), would require 
the knowledge of the orbital separation of these systems. 

It has sometimes been argued (e.g. Eggleton 1986; 
Tout & Eggleton 1988; Kenyon 1994; Han et al. 1995) 
that the mass loss rates of giant stars in binary systems 
must be larger than those of single red giants. More pre- 
cisely, Han et al. (1995) suggest that the mass loss rate 
of a giant approaching its Roche lobe in a binary system 
exceeds by more than a factor 10^ the rate predicted by 
the Reimers formula (Reimers 1975). As far as the binary, 
extrinsic S stars are concerned, this effect, if present, is 
clearly not large enough to bring their mass loss rates to 
the level of the intrinsic S stars. The difference in the evo- 
lutionary stages of extrinsic and intrinsic S stars, believed 
to be RGB (or Early-AGB. according to the terminology 
of Iben & Renzini 1983) and TP-AGB stars, respectively, 
thus appears to be of greater importance for the mass loss 
rate than their binary or non-binary character. 

7. Summary 

Our extensive discussion of the circumstellar properties of 
S stars confirms the dichotomy extrinsic/intrinsic S stars 
and leads to the following conclusions. 

Extrinsic S stars have the lowest mass loss rates among 
S stars (< 2 10~^ Mq y~^), and undetectable circumstel- 
lar shells. This is consistent with the hypothesis that these 
stars are much less evolved than intrinsic S stars, popu- 
lating the RGB or Early-AGB rather than the TP-AGB 
like intrinsic S stars. The binary character of these stars 
does not seem to increase their mass loss rate, as it has 
sometimes been suggested, e.g., in discussions relating to 
symbiotic or barium stars. Note, however, that the few bi- 
nary intrinsic S stars are among the stars with the largest 
mass loss rates in our sample. 

Among intrinsic S stars, various subclasses must be 
distinguished: 

1. A few Tc-rich S stars (e.g., HR Peg) have very low 
mass-loss rates, similar to those of the extrinsic S stars, 
with undetectable circumstellar shells [Region A of the 
{K — [12], [25] — [60]) color-color diagram]. These are S 
stars with weak chemical peculiarities, barely distinguish- 
able from normal M giants. 

2. S stars with weak chemical peculiarities and tenuous 
0-rich circumstellar shells (small K — [12] excess), fed 
by a small albeit measurable mass loss rate of a few 
10~^ Mq (Region B). These are short-period (P ^ 
100 to 150 d) SR variables or short-period (P < 500 d) 
Mira variables. This region may be contaminated by (pos- 
sibly normal M) supergiants of variability type SRc (T 
Get, RS Cnc, Y Lyn). These SRc supergiants have re- 
solved shells at 60 fim, and for RS Cnc, a double wind. 

3. S stars with strong chemical peculiarities and dense O- 

rich circumstellar shells (as indicated by SiO maser emis- 
sion, the 9.7 fim silicate feature and large K — [12] in- 
dices), populating the lower part of Region C. Mass loss 



rates range from several 10 to 10 ^ M© y ^. Almost all 
of them are long-period (P > 300 d) Mira variables. 

4. S stars with strong chemical peculiarities (often classi- 
fied as SG in the optical), with neither SiO maser emission 
nor 9.7 /xm silicate emission, with featureless IRAS spec- 
tra and moderate to large [25] — [60] indices. These stars 
populate Region D and the upper part of Region C (like 
optical carbon stars) and have moderate mass loss rates (a 
few 10~^ to a few 10~^ Mq y~^) with small wind veloci- 
ties (generally ^ 8 km s~^). They are mainly semi-regular 
(with periods ranging from 60 to 360 d) or irregular vari- 
ables, with a few short-period (P < 370 d) Mira variables. 
These properties arc reminiscent of carbon stars. All their 
features may equally well be explained by a carbon-rich 
circumstellar shell, or by a detached shell. The relation- 
ship between stars in Regions C and D is far from being 
clear. It may be a difference of galactic populations, as 
suggested by their different properties as variable stars, or 
stars might oscillate between Regions G and D at differ- 
ent phases of their variability cycle, as observed for some 
variable IRAS sources in our sample, or they may be on 
different parts of the loop associated with the detachment 
of a dust shell. 

5. A few S stars have well-resolved (and thus very ex- 
tended) IR shells (also visible in the optical in the case of 
RZ Sgr), probably detached from their parent star, as in- 
dicated by strong 60 and 100 /^m excesses (Region E). The 
mass loss rates span a wide range, from a few 10~^ to sev- 
eral 10~^ Mq y"-'^, as do the wind velocities (25 km 

for TT Gen to 2.8 km s"^ for FU Mon). These stars with 
extreme properties are often SC or OS stars, and may be 
experiencing a loop in the IRAS color-color diagram as 
first proposed by Willems & de Jong (1988), triggered by 
their C/0 ratio being very close to unity. In TT Gen, the 
dominant spectral features are ZrO bands at some times, 
and G2 bands at other times. This is a rare and short-lived 
phase, and thus does not face the difficulties generally op- 
posed to the Willems & de Jong scenario. Some SG or GS 
stars in Region G, like BH Gru (a sister case of TT Gen) 
and UY Gen, may be at the end of the counter-clockwise 
loop in the IRAS color-color diagram, entering the region 
of heavily-obscured IR carbon stars. Region E is proba- 
bly fed by stars leaving Region C when their C/0 ratio 
reaches values close to unity. However, there is no need 
that all S stars experience such a loop. 
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Table Al. IRAS flux densities for x Cyg 

JD- 2 440 000 F12 offset det F25 F60 [12] -[25] [25] -[60] 
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Appendix A: Stellar variability and the IRAS flux 
densities 

As discussed in Sect. 2.2, the flux densities obtained from 
co-addition of the IRAS data generally agree well with 
those in the PSC with, however, some notable exceptions, 
as shown in Fig. 1. Several stars have flux densities which 
are in disagreement at a level well outside that expected 
from noise. In this Appendix we investigate whether these 
disagreements can be explained by stellar variability, using 
the data for x Cyg as an example. The cphemeridcs listed 
in the GCVS give a variation of about 6 mag at V between 
maximum and minimum light, a period of 408 d and the 
epoch of zero phase which is closest to the IRAS launch 
date of JD 2445404.4. 

The colors of x Cyg, calculated from the IRAS 
PSC flux densities (see Table A2), arc photospheric, de- 
spite the strong circumstellar CO emission (discussed in 
Sect. 5). On the other hand, model IRAS colors calculated 
from the mass loss rate given by the CO lines (several 
10~^ M0 y~^), assuming a normal gas to dust ratio and 
silicate grains (the star is an SiO maser; e.g. Patel et al. 
1992), are well displaced from photospheric values. This 
discrepancy does not mean that the circumstellar envelope 
is dust-free, however; the ratio of the 12/im to 2/um flux 



densities is well in excess of the photospheric value for a 
temperature of 2400 K (Haniff et al. 1995), demonstrat- 
ing the presence of an appropriate amount of circumstellar 
dust. 

We therefore re-examined the raw IRAS data for this 
star. The individual IRAS 12/im, 25/Ltm, and 60/im obser- 
vations for X Cyg are listed in Table Al. Column 1 gives 
the Julian date of the observations computed from the 
'Satellite Operation Plan' number attached to each scan 
and from the mission chronology provided by the Explana- 
tory Supplement. Column 2 lists the phase. Next are the 
flux density, the offset between the scan center and the 
stellar position, and the detector number for the 12 /xm 
data. The final columns contain the 25 /im and 60 /im 
flux densities and the [12] — [25] and [25] — [60] colors cal- 
culated from combinations of the observed flux densities. 
The colors calculated from the flux densities derived from 
individual scans locate x Cyg in Region B of the color- 
color diagram. 

Table Al shows that there are occasionally large dis- 
crepancies between the flux densities mcasiired by differ- 
ent detectors at the same time; see, for example, the 12 /xm 
flux densities observed at JD2445466.1. The low flux den- 
sity is probably due to the large offset between the stel- 
lar and detector center positions; such discrepant observa- 
tions are usually flltered out in the IRAS data processing. 
Real variability is also apparent. The data were obtained 
at two epochs, near maximum and minimum phase, and 
the flux densities at all four wavelengths are systematically 
higher for the first set of observations than the second. The 
mean flux densities for these two epochs are listed in Ta- 
ble A2; the differences are A = (F^ax - -P'min)/-Fmax = 
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32% at 12 nm, 26% at 25 ;um, 25% at 60 /zm and 31% at 
100 Atm. 

Can these variations be attributed to the stellar vari- 
ability? The huge variations in Mira variables at visible 
"wavelengths are due largely to the changing photospheric 
temperature; the variation in the bolometric magnitude 
is much smaller, about a factor of 2 (e.g. Petit 1982; 
Hoffmeister et al. 1985). Because of the variation in stel- 
lar effective temperature, the star reaches maximum light 
later at longer wavelengths in the visible and near infrared 
(e.g. Le Bcrtre 1992). In particular, the variation at visible 
wavelengths leads the total light variation by about 0.1 of 
a period. The two groups of observations in Table Al were 
thus made close to maximum and minimum luminosity. 

We modeled the object as a star with a circumstellar 
envelope and varied the luminosity of the model star. The 
envelope contains silicate grains, has a dust loss rate of 
8 10~^° M0 y~^ and is assumed to have the same out- 
flow speed as the gas (9.5 km s"-'^). The star is assumed 
to be a black body of temperature 2400 K and luminosity 
3000 Lq at minimum, and 2800 K and 6000 at maxi- 
mum. The resulting model flux densities are listed in Table 
A2 and are reasonably close to the observed values (note 
that the discrepancy at 12 /itm may be caused in part by 
the saturation of the detectors at these flux densities ex- 
ceeding 1000 Jy). The variation in the IRAS flux densities 
can thus be fully explained by the stellar variability. 

These results show that caution is required in the inter- 
pretation of IRAS data for variable stars; colors must be 
calculated from data taken at the same epoch, as should 
models of the circumstellar envelope. Variability intro- 
duces a significant amount of scatter into the IRAS colors, 
especially [12] — [25], which as a result does not provide 
the clean measure of the stellar mass loss rate which the 
models predict (see, for example, the discussion by Habing 
1996). 

Table A2. Observed and model IRAS flux densities for x Cyg. 

The PSC fluxes are also listed 

Observed (averages) 



Phase 


F12 


F25 


F60 


FlOO 




(Jy) 


(Jy) 


(Jy) 


(Jy) 


max 


1804 


568 


96 


20 


min 


1223 


420 


72 


14 
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1688 
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81 


18 




Calculated 




Phase 


F12 


F2r, 


F60 


f 100 




(■Jy) 


(Jy) 


(■Jy) 


(Jy) 


max 


I486 


576 


86 


20 


min 


1108 


486 


79 


17 
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Table 1. IRAS co-added flux densities for S stars, grouped according to their location in the {K — [12], [25] — [60]) diagram 
Region A: Stellar photospheres 
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a: Strong cirrus contamination 

Region B: Small ZrO index and no [25] — [60] excess 



GGGSS IRAS 


F2.2 


R.cf 


F12 


F25 


F60 


FlOO 


Te 


LRS 


VG 


Sp 


Var 


P 


AV 


Name 






(•Jy) 




(■Jy) 


(■Jy) 


(.ly) 


(■Jy) 












w 


(mag) 




8 


00192-2020 


1356.0 


2 


199.1 


81.66 


14.57+ 


6.08+ 


yes 


16 




M5-6Se 


SRc 


159 


1.9 


T Get 


134 


05208-0436 


80.2 


8 


20.56 


10.85 


1.65 






43 


U 


M4Swk 








V535 Ori 


168 


05495+1547 


26.1 


1 


11.35 


3.52 


0.62 








F 


S7.5,le 


M 


494 


4.9 


Z Tau 


221 


06266-1148 


27.6 


7 














F 


S-*2e 














,1D 2 445 420: 


8.97 


4.10 


0.94 
























JD 2 445 620: 


14.57 


7.64 


1.14 




















265 


06466-2022 


87.2 


2 


8.34 


2.41 


0.47 






16 




M4S 








HD 49683 


323 


07117-1430 


40.6 


2 


4.17 


1.52 


0.30 




















408 


07461-3705 






14.11 


6.25 


0.87 






14 


F 


S6*l 










436 


07545-4400 


51.6 


1 












22 


E 


S4,2 


M 


340 


4.7 


SU Pup 






JD 2445 470: 


15.44 


6.88 


1.29 
























JD 2 445 645: 


25.03 


11.95 


1.96 




















446 


07573-6509 


27.6 


1 


5.45 


1.72 


0.27 










S7,2 


M 


280 


2.7 


X Vol 


474 


08098-2809 


55.5 


2 


7.07 


2.23 


0.42 








S 


S4,2 








-27:5131 


533 


08348-3617 


55.5 


1 


11.25 


3.44 


0.53 






18 


S 


S5,2 








-36:4827 


589 


09076+3110 


2808.0 


2 


489.88 


210.29 


34.57+ 


11.71+ 


yes 


22 




M6S 


SRc 


120 


1.5 


RS Gnc(OH?) 


626 


09411-1820 






10.0 


4.10 


0.71 








F 


MOS 


M 


300 


5.0 


FM Hya 


704 


10436-3459 


133.2 


1 


31.32 


13.98 


1.93 


1.09 




42 


E 


S5,4 


SR 


104 


2.0 


Z Ant 


903 


15492+4837 


1295.0 


2 


205.41 


98.75 


17.13+ 


7.14+ 


yes 


41 




M6.5S 


SRb 


148 


1.5 


ST Her 


914 


16097-6158 






12.24 


4.51 


0.70 








F 


S4,l 


M 


323 


5.0 


Y TrA 


948 


16552-5335 


39.1 


1 


8.80 


2.90 












S4,4 










1099 


19008+1210 


105.8 


2 


11.54 


3.62 






yes 




S 


S5/2 


Lb 




0.6 


V915 Aql 


1131 


19226-2012 


20.5 


7 


7.32 


2.94 


0.50 








F 


M8Swk 


M 


332 


7.0 


TT Sgr 


1165 


19486+3247 


5813. 


2 










yes 




E 


S7/1.5e 


M 


408 


6. 


X Cyg (SiO) 






JD 2 445 450: 


1781.65" 


552.84 


95.47+ 


21.32 






















JD 2445640: 


1215.31" 


414.40 


71.64+ 


14.50 


















1211 


20213+0047 


152.9 


2 


22.54 


7.27 


1.29 






17 


S 


S7,2 


M 


364 


4.5 


V865 Aql 


1346 


23595-1457 


120.3 


2 


13.11 


3.88 


0.73 




yes 


16 


F 


S5-7/1.5-3e 


M 


351 


7.7 


W Get 



a: detector probably saturated 

+: the 'zero-crossing' flux density Fz has been adopted instead of the template flux density Ft, indicating a possibly resolved 
shell (see Sect. 2.2) 
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Table 1. Continued. 

Region C: [12] — [25] and [25] — [60] excesses, silicate emission common 



GCGSS IRAS F2.2 

(Jy) 


Ref 


F12 

(Jy) 


F25 

(Jy) 


F60 

(Jy) 


FlOO 

(Jy) 


Tc 


LRS 


VC 


Sp 


Var 


P 
(d) 


Ay 

(mag) 


Name 


6 


00135+4644 37.4 


2 


15.96 


6.24 


1.57 


0.85 




16 


F 


S4/7e 


M 


346 


6.9 


X And 


9 


00213+3817 453.3 


2 


335.74 


175.66 


26.41+ 


9.45 


yes 




E 


S5-7/4-5e 


M 


409 


9.1 


R And (SiO) 


1 A 

14 


U044o+3Z/i4 9D.0 


2 


37.75 


18.58 


3.97 


2.10 




22 


F 


S6/2e 


M 


430 


7.8 


RW And 


28 


01159+7220 147.4 


o 
z 


o4o. ( D 


1 no 
iyz.Oz 


Qn 1 A 1 


1U.88 




zz 


T? 

rj 


b4,DC 


iVl 


DlZ 


O.Z 


b Uas (biUj 


ou 




D 


D.oy 




r\ AK 










oD/8e 


iVl 


Q p: p: 
000 


0.0 


KZj l^er 


49 


02143+4404 275.7 


o 
z 


iDD.4z 


CIA 
/ U.y4 


1 A ni 
14. Ui 


0.4o 


yes 


zz 


T? 


k3 < / le 


IVl 


oyb 


/ .y 


W Ana (biwj 


149 


05374+3153 280.8 


2 


45.33 


23.32 


5.00 




yes 


43 




M2S 


Lc 




0.2 


NO Aur 


231 


06331+1415 136.9 


2 


21.95 


10.35 


2.61" 


3.90 


no 


42 


F 


S8,5 


SRa 


1145 


1.4 


DY Gem 


283 


06571+5524 118.1 


2 












16 


F 


S5/5e 


M 


378 


7.1 


R Lyn (SiO) 




JD 2 445 


420: 


17.55 


5.61 


1.25 






















JD 2445 


620: 


28.46 


9.46 


2.01 




















307 


07043+2246 88.0 


2 












16 


F 


S3,9e 


M 


370 


8.0 


R Gem 




JD 2 445 


435: 


22.23 


7.53 


2.21 


1.40 




















JD 2445 


620: 


42.33 


16.81 


3.73 


1.48 


















oiD 


U / Uyz+U i oo 




IZ.ZO 


o.oo 


l.UO 


U.OD 








oe 


ivi 


/1 90 
4ZU 


A ft 
4.0 


VV -\ V^iVll 




nvi /io_i_ni 1 1 vn 


Z 












1 

10 




ivi ( se. 


A/r 
ivi 


oyO 


7 /I 
( .4 


111 V ivion 




T'n o /I /I ^ 

JU Z 440 


44U. 


Z / . ( D 


1 1 7t^ 
11. / 


9 

Z.Ol 


U.oO 




















JD Z 440 


ozo. 


1 Pi 77 


7 91 
/ .Zl 


1 fiQ 


i.UO 


















341 


07197-1451 




20.24 


10.47 


1.86 






27 


E 




M 


314 


4.0 


TT CMa 


347 


07245+4605 964.7 


2 


133.50 


65.92 


13.20+'* 


5.58+ 


yes 


23 




M6S 


SRc 


110 


2.5 


Y Lyn (SiO) 


014 


Uyooo— oo4y Do.U 


1 












01 


E 


S7,8e 


M 


408 


3.0 


UU Vel 






510: 


11.53 


5.98 


1.77 






















-\T\ O AAt^ 
JU Z 440 


DoU. 


1 7 7Q 


y .o ( 






















649 


09564-5837 




156.15 


71.37 


14.31 






15 


E 


S6.5/1- 


SRb 


109 


1.3 


RR Car 


656 


10017-7012 43.7 


1 


6.02 


1.72 


0.37 










S5,6 


M 




3.5 


KN Car 


816 


13136-4426 356.8 


1 


56.55 


20.77 


4.04 


1.78 




43 


C 


S6/8=CS 


SR 


114 


2.0 


UY Can 


ozl 


loZ4:\j—o <4:Z 




13.40 


6.27 


6 






14 


F 


S6*3 


SR 


198 


1.9 


EE Cen 


80 1 


i4o / z— oiUo 4o.o 


9 


21.38 


9.13 








15 


E 










CSV2170 


872 




1 

i 


01 SI 
zl.oi 


Q 

y .zo 


9 18 
Z.lo 








r 


/ ,8e 


AT 
ivi 


OZD 


O.D 


v^i ijup 


931 


16334-3107 377.0 


o 
z 


0z.4y 


Zi.88 


4.ZZ 






10 




b8/4var 


bita 


1 (\A 

iy4 


0.0 


b 1 bCO 


954 


17001-3651 453.3 


1 












22 


E 


S7,2 


M 


449 


8.2 


RT Sm 




JD 2 445 


410: 


170.65 


69.64 


15.37 






















JD 2445 


600: 


269.51 


118.53 


25.01 






















17081+6422 377.0 


z 


Di.yb 


zo.oz 


A 

4.8 ( 


o 07 
A.Li 










IjD 




U.4 


1 V Ura 


1093 


18575-0139 37.7 


X 


17.18 


5.92 


1.83 










SC 


M: 




3.8 


VX Aql 


1096 


18586-1249 242.3 


2 


51.78 


19.19 


4.03 






21 


E 


S6/3e 


M 


395 


8.8 


ST Sgr 


1112 


19111+2555 18.7 


6 












41 




SC 


M 


438 


5.8 


SLyr 




JD 2 445 


440: 


43.81 


21.43 


4.39 






















JD 2445 


630: 


84.17 


38.23 


7.53 


2.73 


















1115 


19126-0708 286.0 


2 












22 


E 


S6/6e 


M 


490 


7.0 


W Aql (SiO) 




JD 2 445 


440: 


1397.52= 


719.87 


132.67+ 


36.87 




















JD 2445 


630: 


1057. 16*^ 


481.93 


100.58+ 


24.90 


















1117 


19133-1703 166.1 


2 


42.27 


14.62 


4.76+ 


4.76+ 


yes 




F 


S5/6e 


M 


394 


5.8 


TSgr 


1150 


19354+5005 107.7 


2 


107.44 


52.86 


12.22+ 


5.52 


yes 


22 


E 


S6/6c 


M 


426 


8.3 


R Cyg (SiO) 


1159 


19451+0827 15.3 


7 


7.69 


3.53 


0.84 








I 


Sc 


M 


607 


2.0 


QU Aql 


1175 


19545-1122 




7.59 


3.80 








29 


E 


M6S 








V1407 Aql 



a: close weak source 

b: profile badly distorted by close source 
c: close weak source apparent on some scans 

d: moves to Region B after removing the extended shell contribution 
e: Detector possibly saturated 

+: the 'zero-crossing' flux density has been adopted instead of the template flux density Ft, indicating a possibly resolved 
sheU (see Sect. 2.2) 
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Table 1. Region C (continued) 

GCGSS IRAS F2.2 Ref F12 F25 F60 FIDO To LRS VC Sp Var P AV Name 

(Jy) (Jy) (Jy) (Jy) (Jy) (d) (mag) 



1200 


20114+7702 












I 


S5/6e 


M 


326 


6.9 


SZ Cep 




JD 2 445 394: 


7.27 


2.71 


0.68 




















JD 2 445 580: 


9.67 


3.45 


0.72 


















1268 


21172-4819 14.9 1 


7.15 


2.34 


0.63 






I 


S2,5 








-48:13866 


C3107" 


22036+3315 44.1 2 


15.92 


7.28 


1.60 






C 


SC 


M 


439 


6 


RZ Peg 


1294 


22196-4612 4211. 1 










yes 42 


E 


S5,7 


SRb 


150 


1.3 


tt' Gru 




JD 2 445 470: 


969.23 


471.63 


88.71+ 


23.88 


















JD 2 445 645: 


1421.10'' 


419.22 


92.25+ 


25.80 
















1345 


23554+5612 66.7 2 


52.07 


28.10 


8.80 


5.28 


42 


E 


S6/6e 


M 


476 


6.9 


WY Gas 


1347 


00001+4826 18.2 6 


50.54 


24.54 


4.15 




21 


E 


S5/6e 


M 


396 


3.3 


IW Gas 



a: Star number from the General Catalogue of Carbon Stars (Stephenson 1989) 

b: Detector possibly saturated 

+: the 'zero-crossing' flux density has been adopted instead of the template flux density Ft, indicating a possibly resolved 
shell (see Sect. 2.2) 
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Table 1. Continued. 

Region D: No silicate emission, weak K — [12] excess, some [25] — [60] excess 



GCGSS IRAS 


F2.2 


Ref F12 


F25 


F60 


FlOO Tc 


LRS 


VC 


Sd 


Var 


p 


AV 


Name 






CTv") 




fJv) 


CTv") 




Civ) 










(d) 






12 


00435+4758 


64.9 


2 


8.47 


2.49 


a 


" VP'S 


01 


s 


S5/3c 


M 


277 


7.7 


U Cas 


17 


00486+3406 


ou.o 


5 


4.56 


1.44 


KJ.OO 










M 




7 9 


RR AnH 


20 


00578+5620 


169.2 


2 


21.87 


7.81 


2.43 




16 


p 


S6/3 


SRb 


136 


1.8 


V365 Ca<) 

V OU<J V_yclo 


29 


01186+6634 


46.6 


2 


5.10 


1.90 


0.42 


















57 


02228+3753 


95.6 


2 


9.89 


2.80 


0.66+ 


0.70 


01 


s 


S8,8 


SR 


159 


1.1 


BI And 


103 


04352+6602 


398.5 


2 


42.56 


11.64 


3.70 


2 21 veq 


17 


s 


S4,7e 


M 


373 


7.1 


T Cam 


116 


04543+4829 


175.5 


2 


12.81 


4.27 


1.38 






F 


S5,8 


SRb 


182 


0.7 


TV Aur 


160 


05440+1753 






4.72 


1.55 


0.50 










SRa 


364 


2.5 


EI Tau 


237 


06347+0057 


51.1 


2 


6.18 


2.21 


0.85 










Lb 




0.9 


CX Men 


312 


070Q5+6S53 


138.2 


2 


14.57 


5.98 


1.94+ 


1.87 yes 






M5S 


Lb 




0.6 


A A r^am 


387 


07399-1045 


1 77 9 


2 


19.48 


7 9^ 


9 ns 




1 s 


CJ 
o 


oo,u 


ortu 




9 8 


O KJ iVKJli 


524 


08308-1748 


65.5 


2 


5.85 


2.00 


a 










M 


250 


1.5 


SZ Pyx 


556 


08461-7051 


53.0 


1 


5.70 


1.49 


0.46 


0.33 






so 


Lb: 




0.4 


UX Vol 


617 


09358-6536 


39.5 


1 


4.00 


1.22 


0.27 








S7,8 










696 


10389-5149 


109.7 


1 


10.61 


3.21 


1.11 




18 


S 


S5,2 


Lb 




0.4 


HP Vel 


787 


12135-5600 






19.78 


6.55 


1.65 




18 




SC6/8e 


M 




4.3 


BH Cru 


803 


12417+6121 


39.5 


3 


4.39 


1.38 


0.35 


0.39 yes 






S3/6e 


M 


226 


5.6 


S UMa 


830 


13440-5306 


129.5 


1 


15.60 


5.66 


2.08 




18 




SC 


Lb 




1.2 


AM Cen 


834 


13477-6009 


387.6 


1 


49.16 


14.57 


3.42 




17 


s 


S8,5 


SR 


307 


3.3 


VX Cen 


904 


15548-7420 


58.1 


1 


6.12 


1.93 


0.65 


0.59 






SC 


SRa 


152 


3.5 


VY Aps 


923 


16209-2808 


71.2 


2 


7.71 


4.10 


1.27 


















940 


16472-6056 


49.2 


1 


5.71 


2.06 


0.72 








SC 


M: 




2.3 


LQ Ara 


978 


17206-2826 


319.4 


2 


39.93 


11.48 


2.72 




17 


s 


S5,4 


SRb 


320 


1.4 


V521 Oph 




17553+4521 


545.0 


2 


55.41 


17.50 


4.14+'^ 


3.20+ yes 


17 




MII-III 


SRb 


120 


0.9 


OP Her 


1011 


17562-1133 






8.45 


3.30 


C 




01 














1056 


18310-3541 


53.0 


7 


5.99 


2.05 


0.62 








S5,4 


L 




0.5 


V3574 Sgr 


1141 


19311+2332 


258.2 


2 


32.48 


13.87 


3.06 




16 


F 


S6/5 


Lb 




3.7 


EP Vul (SiO) 


1188 


20026+3640 


356.8 


2 


40.22 


15.45 


5.37 


yes 


31 


S 


S7,5 


SRb 


212 


3.0 


AA Cyg 


1310 


22479+5923 


91.3 


2 


10.03 


3.25 


a 






F 


S6/2 


SRb 


60 


1.4 


CV Cep 


1308 


22476+4047 


766.3 


2 


100.51 


32.66 


10.25+ 


10.07+ 


16 


S 


S7.5/1 


SRb 


174 


1.7 


RX Lac 


1339 


23380+7009 


60.9 


2 


7.95 


2.89 


0.86 






F 


M6S 








S1339 



a: profile badly distorted by close source 
c: strong cirrus contamination 

d: moves to Region B after removing the extended shell contribution 

+: the 'zero-crossing' flux density has been adopted instead of the template flux density Ft, indicating a possibly resolved 
shell (see Sect. 2.2) 
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Table 1. Continued. 

Region E: Strong [25] — [60] excess, many resolved shells 



GCGSS IRAS 


F2.2 


Ref F12 


F25 


F60 


FlOO 


Tc 


LRS 


VC 


Sp 


Var 


P 


AV 


Name 






(Jy) 




(Jy) 


(Jy) 


(Jy) 


(Jy) 












(d) 


(mag) 




82 


03452+5301 


66.7 


2 


7.66 


3.36 


2.40** 






01 


F 


S5,8 


Lb 




0.8 


WX Cam 


89 


04123+2357 


258.5 


2 


25.13 


7.63 


3.30+ 


2.56 


yes 




S 










+23:654 


114 


04497+1410 


1010. 


2 


87.25 


22.19 


7.52+ 


3.57 


yes 


18 




S3,l 


SRb 


30 


0.2 


0^ Ori 


117 


04599+1514 


88.8 


2 


11.11 


3.45 


1.44 






17 




sc 


SRb 


370 


1.1 


GP Ori 


212 


06197+0327 


144.7 


2 


17.16 


4.63 


2.66+ 


5.88+ 




18 


S 


sc 


SR 


310 


1.1 


FU Men 


422 


07507-1129 


79.5 


2 


7.21 


1.82 


0.48 




yes 




I 


S5/2 


Lb 




0.1 


NQ Pup 


540 


08388-5116 


3.9 


1 


5.39 


2.02 


0.97 










S6,8 








S325 


652 


10015-4634 


58.7 


1 


6.31 


1.84 


1.23 










S5,8 










817 


13163-6031 


94.7 


1 


21.52 


12.91 


7.01 










SC5:/8 








TT Cen 


929 


16316-5026 


71.2 


8 












42 




Se 












JD 


2 445 380: 


65.34 


31.24 


16.70+ 


30.18+ 




















JD 2 445 600: 


98.46 


45.58 


20.72+ 


38.74+ 


















1189 


20044+2417 


129.5 


2 


19.16 


7.55 


4.95+ 


7.30+ 




16 


F 


S4,2 


SRa 


370 


2.0 


DK Vul 


1195 


20100-6225 


133.2 


1 


11.03 


2.82 


1.77+ 


3.27+ 


yes 


18 


S 


S4,4 








HD 191630 


1196 


20120-4433 


187.2 


1 


39.13 


25.58 


13.35+ 


25.88+ 




16 


F 


S4,4 


SRb 


223 


2.6 


RZ Sgr 



b: close weak source apparent on some scans 

+: the 'zero-crossing' flux density has been adopted instead of the template flux density Ft, indicating a possibly resolved 
shell (see Sect. 2.2) 
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Table 2. S stars whose circumstellar envelopes are resolved at 60 /im, listed according to the region to which they belong in 
the {K — [12], [25] — [60]) diagram (sec Fig. 3). In column Rem, YPK+ and YPK— denote stars found to be extended or not, 
respectively, by Young ct al. (1993a). Stars with no YPK label were not examined by Young et al. (1993a). The mass loss rates 
and wind velocities listed in columns M and Vc, respectively, arc taken from Table 6 



IRAS 


Name 


Var 


F60 

(Jy) 


p 


W25 



W50 



QF 


M 

(Ms y-^) 


( km s"^) 


Rem. 














Region B 








00192-2020 


T Get 


SRc 


14.6 


1.05 


1.96 


1.44 


14.5 


4.0(-8)'' 


6.9 


YPK-,b 


19486+3247 


X Cyg 


M 












3.5(-7) 


9.5 


YPK+,f 




JD 


2 445 450: 


95.5 


1.01 


1.96 


1.45 


9.5 










JD 


2 445 640: 


71.6 


1.002 


1.93 


1.42 


0.4 








09076+3110 


RS Cnc 


SRc 


34.6 


1.04 


1.97 


1.44 


29.2 


1.7(-7) 


7.4 


YPK+,b,f 




RS Gnc 














4.5(-8) 


2.6 


e 


15492+4837 


ST Her 


SRb 


17.1 


1.02 


1.95 


1.43 


6.7 


1.7(-7) 


9.5 


YPK+,b 














Rei 


;ion G 








00213+3817 


R And 


M 


26.4 


1.05 


2.01 


1.47 


30.6 


7.5(-7) 


9.3 


YPK- 


01159+7220 


S Gas 


M 


30.1 


1.01 


1.94 


1.42 


10.4 


3.5(-6) 


20.0 


YPK- 


07245+4605 


Y Lyn 


SRc 


13.2 


1.15 


2.00 


1.45 


30.9 


2.0(-7) 


8.6 


YPK+,d 


1 Q1 C\'7r\Q 


W Aql 


M 












l.l(-5) 


ic5. / 


I riS. — 




JD 


2 445 440: 


132.7 


1.006 


1.93 


1.44 


5.8 










JD 


2 445 630: 


100.6 


1.02 


1.95 


1.44 


7.5 








19133-1703 


TSgr 


M 


4.8 


1.11 


2.05 


1.48 


8.0 


3.5(-7) 


14.1 


b 


19354+5005 


RCyg 


M 


12.2 


1.03 


1.96 


1.44 


10.0 


1.5(-6) 


10.5 


YPK+,f 


22196-4612 


TT^ Gru 


SRb 












4.0(-7) 


11.0 


YPK+,f 




JD 


2 445 470: 


88.7 


1.04 


1.98 


1.44 


37.9 










JD 


2 445 645: 


92.2 


1.07 


2.02 


1.49 


38.3 




















Region D 








02228+3753 


BI And 


SR 


0.7 




2.17 


1.58 


2.6 








07095+6853 


AA Gam Lb 


1.9 




2.08 


1.57 


4.1 


2.5(-7) 


17.9 




17553+4521 


OP Her 


SRb 


4.1 


1.15 


2.16 


1.53 


42.4 






b 


22476+4047 


RX Lac 


SRb 


10.2 


1.16 


2.11 


1.52 


23.6 


1.5(-8) 


3.4 


b 














Region E 








04123+2357 


BD+23°654 - 


3.3 


1.20 


2.18 


1.53 


7.2 








04497+1410 


0^ Ori 


SRb 


7.5 


1.48 


2.89 


1.71 


34.8 


< 4.0(-8) 




c 


06197+0327 


FU Mon 


SR 


2.7 


1.40 


2.42 


1.62 


12.6 


9.0(-8) 


2.8 


b 


16316-5026 


S929 






















JD 2 445 380: 


16.7 


1.50 


3.08 


1.81 


24.8 






b 




JD 2 445 600: 


20.7 


1.89 




1.87 


24.1 






b 


20044+2417 


DK Vul 


SRa 


4.9 


1.22 


2.32 


1.60 


25.3 


2.5(-7) 


5.0 


b 


20100-6225 


HD 191630 - 


1.8 


1.35 


2.57 


1.67 


24.0 


< 9.0(-7) 




b,c 


20120-4433 


RZ Sgr 


SRb 


13.3 


1.26 


2.28 


1.58 


62.1 


1.8(-6) 


8.8 


YPK+,b,d 



Notes: 

a: numbers between parentheses refer to power of ten 

b: 100 ^m profile wider than point-source template 

c: listed as extended at 60 nm in the IRAS Small Scale Structure Catalogue 
d: listed as extended at 100 /xm in the IRAS Small Scale Structure Catalogue 
e: two-components wind 

f: circumstellar shell resolved in GO (Sahai 1992; Stanek et al. 1995) 
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Table 3. Masers in S stars 

Star IRAS SiO H2O OH Refs 



Region A 



BD Cam 


03377+6303 


- 


N 


- 


B96 


NZ Gem 


07392+1419 


- 


N 


- 


B96 


HR Peg 


22521+1640 


N 


N 


- 


B96, Z89 


GZ Peg 


23070+0824 


N 


N 


- 


D73, H89 


Region B 












i Uet 


UUlyz-zUzU 


IN 


IN 


AT 

IN 


Aoya, Ayu, u / o, 












xiy4, L//0, tr { i. 


V 000 \JT1 


UozUo-U4oO 




AT 
IN 


AT 
IN 


TVQQ XJQ1 


Zj lau 


uo4yoH-io4y 


IN 




AT 
IN 


jj/o, jyi 


Jrto *^I1C 


uyu/ D+oiiu 


IN 


IN 


V"? 


L)(o, L)(0, Ayz, 












D82, K77, L78, 












L90, P92, R76 


Z Ant 


10436-3459 






N 


H91 


ST Her 


15492+4837 


- 


N 


N 


D73, D76, D82, 












P79 


S Her 


16496+1501 


N 


N 


- 


B77a, B84, B87, 












B90, E88, S81 


X Cyg 


19486+3249 


Y 


N 


N 


A89b, A92, 












B87, B94, C83a, 












C96, D76, J84, 












J85, J87, J91, 












N85, 085, P92 


W Get 


23595-1457 


N 


N 


- 


B77a, B84, L76, 












L78 


Region C 












sx Ana 


UUzio+ooi / 


V" 

Y 


AT 
IN 


AT 

IN 


ID Ma, D * 0, IDOl, 












TiOA RQ(^ r^Qd 

J3y4, jjyD, KjiJo, 












JJ ( 0, r / 0, jyi, 












iy4 


Ti\M A n rl 

ixvv i\na 


UU44D-rOZZ4 


IN 


AT 
IN 


IN 


Jj / /a, Jj(o, L^oOD, 












r^G« T^Vc; T77'Q 












IV M , W ( Z 




UiioyH- / zzU 


v 

I 


AT 
IN 


IN 


rS / /a, r5y4, r ( 0, 












1 y4 


KZj r^cr 


UiZDD+OUoO 


TvT 
IN 






D M a 


w And 


noi /I Q_|_/| /lA/l 


V" 


AT 
IN 


IN 


AoyD, 0/ /a, 00/. 












TiO/1 r^(M\ TlOf^ 

J3y4, oyo, Jjyo, 












T\'7Q T^TT r\Q(\ 

JJ ( 0, IV / / , UoU, 












iy4, W { Z 


NO Anr 






N 


N 


"Rfil RQfi 


R Lyn 


06571+5524 


Y 


N 




B77a, B94, K84, 












S81 


R Gem 


07043+2246 


N 


N 


N 


B77a, B78, B96, 












J91, L78, W72 


RR Mori 


07149+0111 


N 




N 


K77, L78 


TT CMa 


07197-1451 


N 


N 


N 


B90, B96, H91 


YLyn 


07245+4605 


Y 


N 


N 


A90, B75, B81, 












D73, D82 


RR Car 


09564-5837 


N 


N 




H94, L77 


ST Sco 


16334-3107 


N 


N 


N 


B77a, B94, J88, 












L77, W72 


RT Sco 


17001-3651 


N 


N 


N 


B77b, C71, L76, 












L77, L78, S88 



Table 3. Masers in S stars (Continued) 



Star 


IRAS 


SiO 


H2O 


OH 


Refs 




Region 


C (continued) 










TV Dra 


17081+6422 


N 


- 


N 


B94, 


W72 


ST Sgr 


18586-1249 


N 


N 


" 


B77a 




W Aql 


19126-0708 


Y 


N 


N 


B77a, 


B84, B87, 












B94, 


D89, L78, 












W72, 


Z79 


T Sgr 


19133-1703 


N 


N 




B77a, 


B90 


R Cyg 


19354+5005 


Y 


N 


N 


B94, 


B9d, 095, 












P92, 


W72, Z79 


RZ Peg 


22036+3315 


- 


N 




B96 




7T Gru 


22196-4612 


N 


N 




D89, 


H90 


WY Cas 


23554+5612 


- 


N 




B77a, 


C83b 


Region 


1 J 












U Cas 


00435+4758 


N 


N 




B77a, 


C83b, J91, 












S81 




RR And 


00486+3406 


N 


N 


N 


B77a, 


F78 


S29 


01186+6634 


- 


N 




C83b 




BI And 


02228+3753 


- 


N 


- 


C83b 




1 Cam 


04352+6602 


N 


N 




B77a, 


B94, D76, 












J91 




TV Aur 


04543+4829 


N 






B77a, 


B94 


rjL lau 


05440+1753 


N 






-D / / a 




AA Cam 


07095+6853 




N 




B96 




SU Mon 


07399-1045 


N 






B94 




SZ Pyx 


08308-1748 




N 




C83b 




S UMa 


12417+6121 


N 


N 




B77a, 


B96, D78 


VX Cen 


13477-6009 


N 






P92 




OP Her 


17553+4521 




N 




B96 




EP Vul 


19311+2332 


Y 




N 


B94, 


R76 


AA Cyg 


20026+3640 


N 


N 




B77a, 


B94, B96, 












D82 




RX Lac 


22476+4047 




N 


N 


D73, 


D82, W72 


Region 


E 












0^ Ori 


04497+1410 




N 


N 


B96, 


D73, P71 


CP Ori 


04599+1514 


N 






L78 




FU Mon 


06197+0327 


N 






B77a 




RZ Sgr 


20120-4433 






N 


C71 
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References to Table 3. The parentheses at the end of each ref- 
erence indicate which maser is probed by the paper 
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R.M., Wright A. E., 1989, MNRAS 236, 363 (SiO) 
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Table 4. S stars observed in CO at the Caltech Submillimeter Observatory 



GCGSS 


Other 


a(1950) 


5(1950) 




6 


SpT 


Var 


P 
(d) 


Vlsr 
(km s"') 


r.m.s. 

(mK) 


89 


HD 26816 


04 12 23.7 


+23 57 15 


171.6 


-19.0 


S 


Lb : 






17 


117 


GP Ori 


04 59 56.7 


+15 14 58 


186.0 


-15.7 


SC7/8 


SRb 


370 : 


+86, +92 


37 


149 


NO Aur 


05 37 26.8 


+31 53 43 


176.9 


+0.7 


M2IIIS 


Lc 




-7, -3 


13 


212 


FU Mon 


06 19 46.1 


+03 27 00 


206.6 


-4.9 


SC 


SR 


310 


-46, -39 


10 


283 


R Lyn 


06 57 10.8 


+55 24 06 


161.2 


+23.4 


S5/5e 


Mira 


378 


+9, +28 


18 


422 


NQ Pup 


07 50 43.6 


-11 29 41 


230.5 


+8.04 


S5/2 


Lb 






42 


589 


RS Cnc 


09 07 37.8 


+31 10 03 


194.5 


+42.1 


M6IIIaSc 


SRc? 


120 


+7, +11 


23 


626 


FM Hya 


09 41 10.3 


-18 20 43 


252.8 


+25.5 


MOS 


Mira 


300 




17 


704 


Z Ant 


10 43 40.2 


-34 59 17 


275.9 


+20.9 


S5,4 


SR 


104 




10 


796 


HR 4755 


12 27 16.8 


-41 27 35 


298.6 


+21.0 


M2II - III 






-2 


10 


803 


S UMa 


12 41 45.5 


+61 22 00 


124.6 


+56.0 


S3/6e 


Mira 


226 


+5, +18 


10 


816 


UY Gen 


13 13 36.8 


-44 26 27 


307.6 


+17.9 


S6/8,GS 


SR 


114 


-23 


39 



Notes to Table 4: 

1. The spectral type in column 8 is from the GGGSS 

2. Column 10 gives the observed range of optical velocities, taken from the following references: Brown et al. 1990 (GCGSS 149 
and 589), Udry et al. (1997) and Jorissen et al. (1997) (GCGSS 117 and 212), Wilson 1953 (GCGSS 283, 803 and 816), HofHeit 
(1982) (GCGSS 796) 

3. All stars observed in the C0(2-l) line except GP Ori, NQ Pup and UY Gen, which were observed in the CO(3-2) line 
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Table 5. C0(2-l) data for the four S stars detected by the CSO observations 

GCGSS Ico T(poak) V,, 
(K X km s (K) (km s -^) (km s -^) 



212 FUMon 0.63 ± 0.05 0.15 ±0.02 -41.7 ±0.2 2.8 ± 0.3 

283 RLyn 0.89 ±0.19 0.08 ± 0.02 +15.7 ±1.4 9.3 ±1.7 

589 RS One 13.3 ± 0.9 0.78 +7.7 7.4 

1.17 +6.8 2.6 

704 ZAnt 0.73 ±0.09 0.07 ± 0.02 -15.6 ± 1.0 7.4 ±1.3 



Table 6. A compilation of CO data for S stars (columns 3-10; for data prior to 1990, see Loup et al. 1993) and mass loss rates fitting these data (columns 11-! 
assuming = 2500 K, R^, = 2.5 x 10^^ cm, and / = n(CO)/n(H2) = 6.5 x 10"* 







Observations 


Model 


IRAS 


Star 


Line 


HPBW 


rms 




Tmb 




K 


Ref 


D 




dM/dt 


radius 


Line 


HPBW 




^co 








n 


(K) 


it/" 1^™ 
km/sj 


(K) 


(km/sj 


(km/sj 




(pc) 


(km/sj 


/AT — 1\ 

(Me yr j 


(cm) ( ) 




( ) 


(K) 


(K km/sj 


Region A: 

03377+6303 


BD Cam 


1- 

2- 1 


50 

13 


0.02 
0.05 










BL94 
SL95 


92 


(8.5) 


< 1.6 10~* 












12272-4127 


HR 4755 


2—1 


30 


0.01 










this work 


145 


(8.5) 


< 2.3 10~^ 












loo / Z — / lOD 




1-0 


45 


0.01 










SL95 


130 


(8.5) 


< 5.8 10~* 












ZZOZi-t-iD^U 


xlrL reg 


1-0 


50 


0.01 










BL94 


125 


(8.5) 


< 6.3 10"* 












Region B: 

00192-2020 


T Cat 


1-0 


100 


0.14 










K97 






















2-1 


30 


0.05 


:'.f. 


0.28 


+23.9 


h :i 


K97 


130 


6.9 


4.6 10"* 


1.26 10^® 6.5 


2-1 


30 


0.24 


3.32 


*09076+3110 


RS Cnc 


1 — n 


'\'\ 

oo 




7 ^ 
t .o 


n 77 


-r 1 .o 


7 fl 


N92 






















2—1 


ou 




1 ^ ^ 
xo.o 


n 78 

U. I o 


-1-7 7 


7 4 


^~rno '\iiTr\v\r 


9nn 


7 4 

( .4; 


1 q in~^ 

-L .y nj 


9 (SI 1 n^^ 8 7 

^.U ± -LU o. / 


2—1 


ou 


n 7Q 

u. 1 y 


1 n R'^ 

xu.uo 






2—1 


30 






1.17 


+6.8 


2.6 


t n 1 Q ismv K 
Llllo vV \JL JV 




2.6 


5.2 10"* 


1 81 1 0'^'' fi 


2—1 


30 


1.18 


5.61 






3-2 


20 


0.07 


29.8 


4.3 


+6.8 


4.8 


S95 






















3-2 


20 


0.07 


21.9 


2.63 


+6.8 


5.8 


K97 


















09410-1820 


FM Hya 


2—1 




n 09 










"("riic txr^'iT'l/' 
tllib yyKJLtk. 


1 snn 


("8 5^ 


^ 1 fl 1 












10436-3459 


Z Ant 


2—1 


30 


0.01 


0.73 


0.07 


—15.6 


7.4 


1" rl 1 C! TKzriTlf 


900 


7.4 


2.9 10"^ 


3 9fi 10^® 9 4 


2—1 


30 


0.07 


0.84 


15492+4837 


ST Her 


1 — n 


22 




4 Q 


0.24 


o.o 


o.u 




9Qn 




9 n 1 ri"'^ 


9 4^ 10^^ ^ fi 


1 — n 


22 


fl 1 7 

U. X i 


^ 47 






2—1 


1 

-LO 




1 7 4 


-L .UO 




8 8 
o.o 












2—1 


1 ^ 


1.12 


9n 






2—1 


1 ^ 

-LO 


n 1 n 

U. -LU 




U.4:0 


—4.4 


Q 1 
y . -L 
























^—9 




n n? 


^ 1 
o. ± 


n ^4 


—4.4 


12.0 


TCQ7 










^—9 
o ^ 


90 
^u 


0.41 


Q 1 
y . X 


19486+3247 


y CvE 

A, o 


1-0 


23 


0.10 




2.20 


+9.6 


10.2 


H90 


106 


9.5 


2.6 10"''' 


2.85 10^® 18 


1-0 


23 


1.68 


37.79 






1-0 


45 




10.6 




+10 


10 


M90 










1-0 


45 


0.68 


13.45 






1-0 


50 


0.02 


11.8 


0.77 


+9.7 


8.9 


BL94 










1-0 


50 


0.57 


11.14 






1-0 


100 


0.07 


3.6 


0.27 


+9.7 


10.1 


K97 










1-0 


100 


0.16 


3.00 






2-1 


13 


0.24 




6.2 


+10.0 


8.8 


SL95 










2-1 


13 


5.47 


105.2 






2-1 


25 


0.03 


25.4 


1.77 


+9.9 


9.0 


BL94 










2-1 


25 


2.45 


43.76 






2-1 


30 


0.15 


22.7 


1.63 


+9.9 


9.9 


K97 










2-1 


30 


1.85 


32.59 






3-2 


20 


0.19 


63.0 


4.7 


+7.8 


9.9 


Y95 










3-2 


20 


2.81 


49.28 






3-2 


20 


0.19 


41.5 


3.0 


+10.6 


9.2 


S95 






















3-2 


20 


0.31 


44.3 


3.37 


+10.4 


9.8 


K97 






















4-3 


15 


0.24 


52.7 


4.0 


+9.9 


9.0 


Y95 










4-3 


15 


3.55 


63.02 



Table 6 (continued). 







Observations 


Model 


IRAS 


Star 


Line 


HPBW 


rms 


^co 


Tmb 


Vc 




Ref 


D 


Ve 


AM/dt 


radius 


Line 


HPBW 


Tmb 


^co 








// 




(K IfTTl /si 
I IV Kill/ o I 




I Kill / o I 


(Kin / b I 






firm /a^ 


/■1V/T„ irr~^1 

(,ivl0 yr ) 


(cm.) { ) 




\ ) 




(K Irm /e\ 


Region B (continued): 
































20213+0047 


V8o5 Aql 


2—1 


25 


0.03 










h5L94 


840 


(8.5) 


< 3.0 10 












23595-1457 


W Get 


1-0 


23 


0.19 


_ 


— 


— 


_ 


H90 


950 


(8.5) 


< 3.4 lO"'' 
















1—0 


50 


0.02 


— 


— 


— 


— 


BL94 






















1-0 


45 


0.04 


— 


— 


— 


— 


SL95 






















3—2 


20 


0.03 


- 


- 


- 


- 


Y95 


















Region C: 




































00213+3817 


R And 


1-0 


45 


— 


3.71 


— 


-16 


11 


M90 


490 


9.3 


8.6 lO"'^ 


5.66 10^*^ 7.7 


1-0 


45 


0.22 


3.71 






1-0 


50 


0.04 


5.73 


0.44 


-15.8 


9.3 


BL94 










1-0 


50 


0.18 


3.02 






2—1 


25 


0.05 


6.5 


0.73 


—15.6 


8.5 


BL94 










2—1 


25 


0.85 


11.46 


01159+7220 


S Gas 


1-0 


23 


0.36 


- 


0.44 


-30.8 


22.0 


H90 


860 


20.0 


4.0 10"^ 


1.0 10^^ 8.0 


1-0 


23 


0.43 


15.6 






1—0 


45 




3.43 




—29 


18 


M90 










1-0 


45 


0.12 


4.29 






1-0 


33 


0.11 


6.0 


0.23 


-26.0 


19.0 


N92 










1-0 


33 


0.22 


7.8 






1 n 
i — u 


50 


0.03 


4.63 


0.14 


-31.0 


22.0 


BL94 










1-0 


50 


0.10 


3.49 






2-1 


25 


0.03 


8.1 


0.29 


—30.0 


19.7 


BL94 










2—1 


25 


0.35 


11.62 


02143+4404 


W And 


1-0 


33 


0.09 


4.5 


0.43 


-35.0 


10.4 


N92 


630 


8.5 


8.0 lO"'' 


5.65 10^® 5.9 


1-0 


33 


0.27 


4.15 






1-0 


50 


0.02 


1.37 


0.12 


-34.2 


7.8 


BL94 










1-0 


50 


0.12 


1.83 






2-1 


25 


0.01 


3.0 


0.29 


-34.4 


6.5 


BL94 










2-1 


25 


0.52 


7.04 






3-2 


20 


0.08 


7.0 


0.59 


-36.7 


8.9 


¥95 










3-2 


20 


0.55 


7.43 


05374+3153 


NO Aur 


2-1 


13 


0.05 










SL95 


620 


(8.5) 


< 1.1 10-'' 
















2-1 


30 


0.01 










this work 


















06331+1415 


DY Gem 


1-0 


50 


0.04 










BL94 


890 


8.0 


2.4 10"'^ 


2.87 10^® 2.1 


1-0 


50 


0.01 


0.18 






2-1 


13 


0.05 




0.25 


-16.0 


8.0 


SL95 










2-1 


13 


0.25 


3.50 


06571+5524 


R Lyn 


1-0 


50 


0.01 


0.32 


0.03 


+15.3 


7.5 


BL94 


950 


8.0 


3.5 10"'^ 


3.51 10^** 2.5 


1-0 


50 


0.017 


0.27 






2-1 


25 


0.01 


0.54 


0.04 


+16.6 


7.7 


BL94 










2-1 


25 


0.092 


1.25 






2-1 


30 


0.02 


0.94 


0.08 


+15.7 


9.3 


this work 










2-1 


30 


0.064 


0.87 



Table 6 (continued). 







Observations 


Model 


IRAS 


Star 


Line 


HPBW 


rms 




Tmb 




K 


Ref 


D 




dM/dt 


radius 


Line 


HPBW 


Tmb 


Ico 








II 


(K) 


(K km/s) 


(K) 


(km/s) 


(km/s) 




(pc) 


(km/s) 


(Mo yr-i) 


(cm) 


(") 






(") 


(K) 


(K km/s) 


Region C (continued): 




































07043+2246 


R Gem 


1—0 


23 


0.10 




0.65 


-59.2 


6.0 


H90 


850 


5.4 


2.3 10"'^ 


7.63 10^^ 


60.0 


1- 





23 


0.10 


1.07 






2-1 


13 


0.12 


- 


1.40 


-59.1 


4.8 


SL95 












z— 


1 
1 


13 


0.55 


4.99 






3-2 


20 


0.07 


1.7 


0.28 


—59.0 


5.5 


K97 












Q 

o— 


o 
z 


20 


0.18 


1.57 


07149+0111 


RR Mon 


1-0 


45 


0.01 










SL95 


1200 


(8.5) 


< 7.0 lO"'' 
















07245+4605 


YLyn 


1-0 


33 


0.08 


2.2 


0.20 


—0.7 


8.6 


N92 


330 


8.6 


2.3 10~^ 


2.73 10^® 


5.5 


1- 





33 


0.14 


2.50 


09338—5349 


UU Vel 


1-0 


45 


0.11 


- 


- 


- 


- 


SL95 


1260 


(8.5) 


< 4.0 10"® 
















13136-4426 


UY Gen 


2-1 


23 


0.02 




0.07 


—28.6 


13.1 


SL95 


550 


13.1 


2.1 10 


Z.ob iU 


o n 

z.y 


o 

z— 


1 

i 


23 


0.07 


1.67 






3—2 


20 


0.04 










this work 






















16334-3107 


ST Sco 


1-0 


50 


0.02 


0.61 


0.06 


-4.5 


6.5 


BL94 


540 


7.1 


1.2 10"'^ 


l.yo iU 


Z.O 


1 

i — 


U 


50 


0.01 


0.15 






1-0 


45 


0.02 




0.14 


— 1.2 


7.3 


SL95 












1- 





45 


0.01 


0.19 






2-1 


25 


0.02 


1.17 


0.14 


-5.8 


7.6 


BL94 












2- 


1 


25 


0.09 


1.24 


17001-3651 


RT Sco 


1-0 


45 


0.02 




0.20 


—44.5 


13.1 


SL95 


490 


11.0 


4.6 lO"'^ 


o. / iU 


O.U 


1 — 


U 


45 


0.06 


1.36 






2-1 


23 


0.05 




0.18 


—AAA 


9.6 


SL95 












O 

z— 


1 


23 


0.31 


5.92 


17081+6422 


TV Dra 


1-0 


50 


0.02 


0.13 


0.03 


+22.4 


3.7 


BL94 


540 


5.0 


4.6 10"® 


1.34 10^® 


1.6 


1- 





50 


0.003 


0.034 






2-1 


13 


0.04 


1.0 


0.15 


+22 


6 


093 












2- 


I 


13 


0.17 


1.67 






2-1 


zo 


U.Ul 


n /in 
U.4U 


U.Uo 


_l01 1 
-\-Zl.L 


o.u 














2- 


1 


zO 


U.Uo 


U.4D 






2-1 


13 


0.09 




0.30 


+21.2 


5.1 


SL95 






















18575-0139 


VX Aql 


1-0 


45 


0.12 










SL95 


1400 


7.8 


3.2 10"'' 


3.41 10^® 


1.6 


1- 





45 


0.01 


0.15 






2-1 


23 


0.01 




0.05 


+6.7 


7.8 


SL95 












2- 


1 


23 


0.05 


0.67 


18586-1249 


ST Sgr 


1-0 


50 


0.01 


0.42 


0.03 


+54.6 


9.3 


BL94 


670 


10.0 


3.5 10"'' 


3.25 10^® 


3.2 


1- 





50 


0.022 


0.04 






1-0 


45 


0.01 




0.06 


+55.1 


9.6 


SL95 












1- 





45 


0.027 


0.54 






2-1 


25 


0.04 










BL94 












2- 


1 


25 


0.12 


2.06 






2-1 


23 


0.03 




0.11 


+49.8 


11.2 


SL95 












2- 


1 


23 


0.14 


2.53 


19111+2555 


SLyr 


2-1 


13 


0.04 




0.40 


+51.2 


13.9 


SL95 


1500 


13.9 


2.1 10"® 


8.10 10^® 


3.6 


2- 


1 


13 


0.40 


9.0 



Table 6 (continued). 





Observations 


Model 


IRAS Star 


Line 


HPBW 


rms 




7m B 




K 


Ref 


D 




dM/dt 


radius 


Line 


HPBW 


Tmb 








II 


(K) 


(K km/s) 


(K) 


(km/s) 


(km/s) 




(pc) 


(km/s) 


(Mq yr-^) 


(cm) 


I ) 




(") 


(K) 


(K km/s) 


Region C (continued): 


































19126—0708 W Aql 


1-0 


23 


0.09 




2.31 


-24.6 


19.3 


H90 


610 


18.7 


1.3 10"^ 


2.12 lO^'' 


23.2 


1-0 


23 


2.98 


102.4 




1-0 


45 


- 


26.2 


- 


-25 


20 


M90 






















1-0 


45 


0.05 


28.4 


0.96 


—24.0 


18.1 


N92 












1-0 


45 


1.14 


35.54 




1-0 


50 


0.03 


27.7 


0.98 


—24.0 


17.7 


BL94 












1—0 


50 


0.95 


29.5 




1-0 


45 


0.22 




1.30 


-24.1 


18.3 


SL95 






















2-1 


25 


0.04 


43.1 


1.60 


-24.1 


17.4 


BL94 












2-1 


25 


26.7 


78.2 




2-1 


30 


0.03 


54.8 


2.10 


-25.2 


19.3 


K97 












2-1 


30 


1.96 


56.98 




2—1 


13 


0.16 


- 


6.05 


-24.0 


18.8 


SL95 












2-1 


13 


6.98 


217.3 




2-1 


23 


0.25 




2.04 


—24.3 


18.0 


SL95 












2—1 


23 


3.05 


90.0 




3-2 


20 


0.10 


93.2 


3.51 


—23.8 


19.6 


K97 












3—2 


20 


2.70 


76.95 


19133-1703 T Sgr 


1-0 


50 


0.02 










BL94 






















1-0 


45 


0.02 


- 


- 


- 


- 


SL95 






















2-1 


zo 


U.U4 










JDJ_iy4 






















2-1 


23 


0.02 




0.06 


+10.8 


14.1 


SL95 


810 


14.1 


4.0 10"'^ 


3.22 10^® 


2.7 


2-1 


23 


0.06 


1.53 


19354+5005 R Cyg 


1-0 


50 


0.01 


2.52 


0.15 


-17.5 


9.9 


BL94 


900 


10.5 


1.7 10"*^ 


8.16 lO" 


6.1 


1-0 


50 


0.11 


2.0 




2-1 


25 


0.02 


5.04 


0.33 


-17.2 


11.4 


BL94 












2-1 


25 


0.43 


7.06 




3-2 


20 


0.17 


14.6 


1.3 


-18.2 


10.4 


S95 












3-2 


20 


0.44 


7.18 


22036+3315 RZ Peg 


2-1 


13 


0.04 




0.21 


-23.4 


12.6 


SL95 


1560 


12.6 


1.1 10"^ 


5.93 10^® 


2.5 


2-1 


13 


0.21 


4.57 


*22196-4612 tt^ Gru 


1-0 


45 


0.05 


6.9 


0.33 


-13.2 


19.4 


N92 


160 


11.0 


4.6 10"'^ 


3.70 10" 


15.5 


1-0 


45 


0.53 


11.80 




1-0 


45 


0.08 




0.60 






SL95 






















2-1 


23 






2.20 


-12.5 


11.0 


S92 












2-1 


23 


2.20 


45.67 




2-1 


23 


0.06 




2.00 






SL95 




















23554+5612 WY Gas 


2-1 


13 


0.04 




0.50 


+7.5 


13.3 


SL95 


1300 


13.3 


1.8 10"® 


7.65 10" 


3.9 


2-1 


13 


0.50 


1.08 



Table 6 (continued). 







Observations 


Model 


TP 




Line 


HPBW 


rms 


^co 


TiviB 


Vc 




Ref 


D 




dM/dt 


radius 


Line 


HPBW 


Tmb 


^co 










(K) 


(K km/s) 


(K) 


(km/s) 


(km/s) 




(pc) 


(km/s) 


(Mo yr-i) 


(cm) (") 




(") 


(K) 


(K km/s) 






































00578+5620 


V365 Cas 


1-0 


50 


0.01 










BL94 






















2-1 


25 


0.03 


- 


- 


- 


- 


BL94 






















2-1 


13 


0.02 


— 


0.07 


-2.1 


7.2 


SL95 


800 


7.2 


7.0 10"® 


1.53 10^® 1.3 


2-1 


13 


0.07 


0.98 


U4ooz-|-DDUz 


1 i^am 


1-0 


50 


0.02 










BL94 






















2-1 


25 


0.05 


- 




- 


- 


BL94 


520 


(8.5) 


< 2.2 10"^ 












04543+4829 


TV Aur 


1-0 


50 


0.02 


- 




- 


- 


BL94 






















2—1 


25 


0.03 










BL94 


780 


(8.5) 


< 2.6 lO"'^ 












07095+6853 


AA Cam 


2-1 


13 


0.02 


- 




-43.5 


17.9 


SL95 


880 


17.9 


2.9 10"^ 


2.67 10^*^ 2.0 


2-1 


13 


0.06 


2.14 


07399-1045 


SU Mon 


1-0 


50 


0.01 


- 


- 


- 


- 


BL94 


780 


(8.5) 


< 2.6 10"'^ 
















1-0 


45 


0.01 










SL95 






















2-1 


25 


0.03 










BL94 


















12417+6121 


S UMa 


2—1 


30 


0.01 










this work 


1650 


(8.5) 


< 4.8 10"''' 












*i Q/1 /in c:Qn(^ 




1-0 


45 


0.02 




0.08 


—27.4 


5.4 


SL95 


910 


5.4 


3.5 10" 


4.13 10 3.0 


1-0 


45 


0.054 


0.55 






2-1 


23 


0.04 




0.20 






SL95 










2-1 


23 


0.27 


2.33 


13477-6009 


VX Cen 


1-0 


45 


0.03 










SL95 


530 


(8.5) 


< 4.6 10"^ 












17206-2826 


V521 Oph 


1-0 


50 


0.01 


- 


- 


- 


- 


BL94 


580 


(8.5) 


< 2.5 10"''' 
















1-0 


45 


0.01 










SL95 




















EP Vul 


1-0 


50 


0.02 


0.94 


0.12 


-0.2 


5.2 


BL94 


650 


5.4 


2.9 10"''' 


3.70 10^® 3.8 


1-0 


50 


0.07 


0.68 






1-0 


45 


0.01 




0.08 


-0.7 


5.6 


SL95 










1-0 


45 


0.08 


0.84 






2-1 


25 


0.02 


1.56 


0.20 


+0.1 


5.5 


BL94 










2-1 


25 


0.34 


2.99 






2-1 


13 


0.14 




0.52 


-0.1 


4.7 


SL95 


















20026+3640 


AA Cyg 


1-0 


50 


0.03 


1.51 


0.11 


+22.7 


9.4 


BL94 


550 


7.1 


3.1 10"'' 


3.46 10^® 4.2 


1-0 


50 


0.06 


0.80 






2-1 


13 


0.13 




0.95 


+27.9 


4.8 


SL95 










2-1 


13 


1.05 


12.82 


22476+4047 


RX Lac 


1-0 


45 


0.13 










SL95 






















2-1 


13 


0.12 




0.68 


-15.4 


3.4 


SL95 


380 


3.4 


4.0 10"® 


1.4 lO" 4.1 


2-1 


13 


0.67 


4.42 



Table 6 (continued). 







Observations 


Model 


IRAS 


Star 


Line 


HPBW 


rms 


I CO 


Tmb 


K 




Ref 


D 


v; 


dM/At 
(Mo yr-i) 


radius 


Line 


HPBW 


Tmb 


^co 










(K) 


(K km/s) 


(K) 


(km/s) 


(km/s) 




(pc) 


(km/s) 


(cm) 


(") 




(") 


(K) 


(K km/s) 


Region E: 






































04123+2357 


GCGSS89 


2-1 


30 


0.02 










this work 


650 


(8.5) 


< 1.7 lO"'^ 














04497+1410 


0^ Ori 


1-0 


45 


0.02 


- 


- 


- 


- 


SL95 


330 


(8.5) 


< 2.0 10"'^ 














04599+1514 


GP Ori 


3-2 


20 


0.04 










this work 


1100 


(8.5) 


< 5.2 lO"'' 














*06197+0327 


FU Mon 


2-1 


30 


0.01 


0.66 


0.16 


-41.7 


2.8 


this work 


860 


2.8 


1.0 10"'^ 


7.70 10" 


59.8 
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References: 

H90: Heske 1990; M90: Marguhs et al 1990; S92: Sahai 1992; N92: Nyman et aL 1992; N93: Nyman et al. 1993; 093: Omont et al. 1993; KJ94: Kahane & Jura 1994; 
BL94: Bieging & Latter 1994; SL95: Sahai & Liechti 1995; Y95: Young 1995; S95: Stanek et al. 1995; K97: Knapp et al. 1997b 

Notes: 

1. 06197+0327, FU Mon: Data from Sahai & Liechti (1995) not used - see text; 

2. 09076+3110, RS Cnc: C0(2-l) hnc profile has two components see text; 

3. 13440-5306, AM Gen: Values for central velocity Vc and outflow velocity Vi given by SL95 are averages of the values for the CO(l-O) and CO(2-l) lines; 

4. 22196-4612, tt^ Gru: SL95 do not list values for Vc and because of the line wings 
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